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BOUSSINESQ/BOUSSINESQ SYSTEMS FOR INTERNAL WAVES WITH A 
FREE SURFACE, AND THE KDV APPROXIMATION 

Vincent Duchene 1 

Abstract. We study here some asymptotic models for the propagation of internal and surface waves in 

'ft ' a two-fluid system. We focus on the so-called long wave regime for one-dimensional waves, and consider 

the case of a flat bottom. Following the method presented in [6] for the one-layer case, we introduce a 
new family of symmetric hyperbolic models, that are equivalent to the classical Boussinesq/Boussinesq 
system displayed in [11]. We study the well-posedness of such systems, and the asymptotic convergence 
of their solutions towards solutions of the full Euler system. Then, we provide a rigorous justification 
of the so-called KdV approximation, stating that any bounded solution of the full Euler system can 
be decomposed into four propagating waves, each of them being well approximated by the solutions 
of uncoupled Korteweg-de Vries equations. Our method also applies for models with the rigid lid 

pvi , assumption, using the Boussinesq/Boussinesq models introduced in [7] . Our explicit and simultaneous 

decomposition allows to study in details the behavior of the flow depending on the depth and density 

yp ' ratios, for both the rigid lid and free surface configurations. In particular, we consider the influence of 

the rigid lid assumption on the evolution of the interface, and specify its domain of validity. Finally, 
solutions of the Boussinesq/Boussinesq systems and the KdV approximation are numerically computed, 

C*) ' using a Crank-Nicholson scheme with a predictive step inspired from [3,4]. 
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1. Introduction 

1.1. Motivation of the problem 

This paper deals with different asymptotic models for the propagation of weakly nonlinear internal waves in a 
two-fluid system. The system we study consists in two layers of immiscible, homogeneous, ideal, incompressible 
and irrotationnal fluids under the only influence of gravity. Since we are interested in KdV equations, which 
are unidirectional, we focus on the one-dimensional case, and the bottom is assumed to be flat (see Figure 1). 

Let us denote by p\ the density of the upper fluid and d\ its depth, pi the density of the lower fluid, and 
di its depth, a\ the typical amplitude of the deformation of the surface and a-i the one of the interface, and 
finally A is a characteristic horizontal length, say the typical wavelength of the interface. Then the regime of 



Keywords and phrases: Internal waves, free surface, rigid lid configuration, long waves, Korteweg-de Vries approximation, 
Boussinesq models. 

1 Departement de Mathematiques et Applications, UMR 8553, Ecole normale superieure, 45 rue d'Ulm, F 75230 Paris cedex 
05, France ; e-mail: duchcnc@dma.ens.fr 

© EDP Sciences, SMAI 1999 



TITLE WILL BE SET BY THE PUBLISHER 




Figure 1. Sketch of the domain. 



the system is described by the following dimensionless parameters: 
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The governing equations of such a system, that we call "full Euler", have been obtained in [16]; we briefly recall 
the system in Section 1.3 below. This system is strongly nonlinear, and remains complex, for its direct study as 
well as for numerical computations. That is why the derivation of approximate asymptotic models has attracted 
lots of interests in the past decades. We focus here on the so-called long wave regime, where the dimensionless 
parameters ei, 62, and fx are small and of the same order: 
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The long wave regime for the one-fluid system (water wave) has been considerably studied, and has led to 
many approximate equations. Among them, of particular interest are the Boussinesq systems (from the work of 
Boussinesq [8,9]), and the KdV approximation (from Korteweg and de Vries [27]). The latter model states that 
any solution of the one-layer water wave problem in the long- wave limit splits up into two counter-propagating 
waves, each of them evolving independently as a solution of a KdV equation. A justification of such models 
has been investigated among others by Craig [13], Schneider- Wayne [38] and Bona-Colin-Lannes [6]. The study 
of internal waves followed quickly. When the surface is assumed to be fixed as a rigid lid, the KdV equations 
related to such a system have been formally obtained by Keulegan [25] and Long [31]. The related Boussinesq- 
typc systems have been justified (among many other asymptotic models) by Bona-Lannes-Saut [7]. When the 
surface is not rigid and allowed to move as a free surface, it is known that there exist two different modes of wave 
motion, corresponding to different linear phase speeds (sec Kakutani-Yamasaki [22], Leone-Segur-Hammack [30], 
Michallct-Barthclcmy [34] and Craig-Guyenne-Kalisch [14] for example). Accordingly, the KdV approximation 
states that any deformation of the surface and/or the interface will split up into four waves, each of them being 
lead by KdV equations. 
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However, to our knowledge, the four KdV equations related to the problem have never been revealed simul- 
taneously, neither rigorously justified with a convergence theorem. Such a precise and complete decomposition 
allows to directly compare the different models. In particular, as an application of our results, we present an 
in-depth study of the influence of the rigid lid assumption on the evolution of the interface, and therefore assert 
the domain of validity of such a hypothesis. 

An intermediate result for the construction of the KdV approximation is the full justification (in the sense of 
well-posedness, and convergence of the solutions towards solutions of the full Euler system) of symmetric coupled 
models, that are equivalent to the Boussinesq-type models derived by Camassa-Choi [11], Craig-Guyenne- 
Kalisch [14] and Duchene [16]. The latter systems are justified only by a consistency result, and the stronger 
properties of the symmetric models make them interesting by their own. 

The construction and rigorous justification of the symmetric Boussinesq-type models and the KdV approxi- 
mation are the main motivation of this article. 

1.2. Main results and outline of the paper 

Our study applies the methods of [6] to the case of internal waves, and accordingly uses as direct supports 
the full Euler system and the Boussinesq/Boussinesq model obtained in [16]. The derivation of the full Euler 
system as governing equations of our problem is quickly explained in Section 1.3. The full Euler system is 
consistent with the Boussinesq/Boussinesq model at order C(/i 2 ), provided that e 1: e 2 — O(^); this result is 
recalled and precisely stated in Proposition 2.1, below. 

The first step of our study lies in the construction of symmetric systems, obtained from the original Boussi- 
nesq/Boussinesq model by using a first order symmctrizcr, and withdrawing the 0(e 2 ) terms. Section 2 is 
dedicated to the construction, analysis and justification of such models. The systems we obtain are of the form 

S + e(Si(U) - S 2 dl])d t U + (e + e(Ei(tf) - ^ 2 d 2 x ))d x U = 0, 

with U(t,x) £ K , symmetric matrices So, So, S*2, £2 € -M^M.) (So and S2 being definite positive), and linear 
mappings S\(-) and Ei(-) with values in real symmetric matrices. The well-posedness of such systems over 
times of order 0(1/ e) is stated in Proposition 2.6. The convergence of their solutions towards solutions of the 
full Euler system is then proved to be of order 0(e), for times of order 0(1/ e), in Proposition 2.8. 

From these models, using the classical WKB method, we prove in Section 3 that a rougher approximation 
consists in four uncoupled KdV equations, that is to say that any bounded solution of the full Euler system in 
the long-wave limit splits up into four propagating waves, two of them moving to the right at different speed, 
and the other two moving to the left, each of them approximated by independent solutions of KdV equations. 
Our main result is Theorem 3.1, which states explicitly the decomposition, and precise the convergence rate 
between bounded solutions of the full Euler system and the solutions of the KdV approximation. The case of 
the rigid-lid configuration can be treated in the same way, and is quickly tackled in Section 3.3. 

The complete, simultaneous decomposition of the flow, and its rigorous justification, is a key point for the 
comparison with different models. In particular, we present in Section 3.4 an in-depth study of the behavior of 
the flow predicted by the KdV approximation, depending of the density ratio 7 and the depth ratio S, for both 
the rigid lid and free surface configurations. As a result, we show that the rigid lid hypothesis is valid only 
for small density differences between the two fluids. To our knowledge, this fact has never been established, 
thought it has been addressed for example in [25,34]. 

Finally, the Boussinesq/Boussinesq models and the KdV approximation are numerically computed and com- 
pared for different initial data and parameters in Section 4. The numerical schemes we use are based on 
Crank-Nicholson methods, with a predictive step in order to deal with the nonlinearities. Such a scheme has 
been introduced by Besse and Bruneau in [3,4]; it is formally of order 2 in space and time, and appears to be 
unconditionally stable. The precise schemes in our framework are presented in Section 4.1. 
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1.3. The full Euler system 

Let us recall here briefly the governing equations of our system (see [16] for more details). The velocity 
fields of an irrotational flow can be expressed as gradients of potentials (that we call <pi for the upper fluid 
and <j>2 for the lower fluid). The conservation of mass, together with the incompressibility of the fluids, yields 
Laplace equation for the potentials. The momentum equations of the Euler equations can then be integrated, 
which yields the Bernoulli equation in terms of potentials. The system is finally closed by kinematic boundary 
conditions (stating that no particles of fluid cross the bottom, the surface or the interface) and pressure laws 
(the pressure is assumed to be constant at the surface, and continuous at the interface, see Remark 1.2 below). 
Thanks to an appropriate scaling, the two-layer full Euler system can be written in dimcnsionless form. 

The key point is then to remark that the system can be reduced into four evolution equations coupling 
Zakharov's canonical variables [42] , specifically the deviation of the free surface and interface from their rest 
position (respectively (i, ( 2 ), and the trace of the velocity potentials of the upper fluid evaluated at the surface 
(namely ip\) and of the lower fluid evaluated at the interface (namely ip2 )• The system is then given by 

ad t Ci £1(^1,^2) = 0, 

1 M 
dt( 2 - -G 2 ^ 2 = 0, 

H £2 (!) 

dtd x ipi + a<9 x £i + —d x {\d x ^i\ 2 ) - /J,e 2 d x Afi = 0, 

d t (d x ^ 2 - 7#(V>i,^2)) + (1 - J)d x ( 2 + jd x (\d x ,p 2 \ 2 



1 \H{i> 1 ,il> 2 )\ 2 )-^ 2 d x N 2 = Q, 



with a = — and where A/i and Af 2 are given by the formulae 



2(l + M |eACi| 2 ) : 



A/' 2 



(iG 2 -0 2 + e2d x ( 2 d x tp2) 



7(7^2 + eid&Htyu V> 2 )) 2 



2(1 + A /|e 2 6>^C2| 2 ) 



and with G\ and G 2 the Dirichlet-Neumann operators and H the interface operator, defined as follows. 
Definition 1.1. For Ci, ( 2 £ W 2 -°°(R) and d x ip!, d x ip 2 £ H 1 / 2 (R), the operators G lt G 2 and H are defined by 

G 2 [e 2 C2]V'2 = (d z <j> 2 - ne 2 d x ( 2 d x <f> 2 )\ z=t2(2 . 
Gi[eiCi, e 2 C2]W'i. ^2) = (dz<f>i - ^id x C\d x 4>i)\ z=l+iiCi , 



with 4>i and 



!>2 the unique solutions of the boundary problems 
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Here, we denote by 

Sl\ = {(x,z) £ M. 2 ,e 2 ( 2 (t,x) < z < 1 + eiCi (*,»)} and 2 = {(x,z)e 

lid 2 + d 2 the scaled Laplace operator, and by d n the upward conormal 
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The domains of the two fluids are assumed to remain strictly connected, i.e. there exists h m [ n > such that 

Vi€l, hi(x) = l + ei77i(a;) > h min > and h 2 {x) = - + e 2 n 2 (x) > h min > 0. (3) 

o 

This assumption is necessary in order to obtain the consistency of the full Euler system (1) with the Boussi- 
nesq/Boussinesq model (4), as seen in [16], and recalled in Proposition 2.1 below. We do not always precise 
this assumption thereafter. 

Remark 1.2. Even if the Cauchy problem associated to the Euler system at the interface of two fluids of 
different positive densities is known to be ill-posed in Sobolev spaces in the absence of surface tension (as 
Kelvin- Hclmholtz instabilities appear), Lannes [28] proved thanks to a stability criterion that adding a small 
amount of surface tension guarantees the well-poscdness of such a problem, with a time of existence that does 
not vanish as the surface tension goes to zero, and thus is consistent with the observations. The stability 
criterion states that the Kelvin-Hclmholtz instabilities appear for high frequencies, where the regularization 
effect of the surface tension is relevant, while the main profile of the wave that we want to capture is located in 
lower frequencies, and is unaffected by surface tension. Consequently, we decide to neglect the surface tension 
term, as its effect does not appear in our asymptotic models. 

Furthermore, we know from Theorem 5.8 of [28] that in the long wave regime, there exists uniformly bounded 
solutions of the full Euler system for times of order 0(1/ fi) (again with a small amount of surface tension) in 
the rigid lid case. With this result in mind, we assume in the following that smooth, uniformly bounded families 
of solutions to (1), whose existence is assumed in Proposition 2.8 and in Theorem 3.1, indeed exist. 

2. Derivation and analysis of Boussinesq/Boussinesq models 

As said in the introduction, the starting point of our study is the Boussinesq/Boussinesq model, obtained 
from (1) thanks to an asymptotic expansion of the operators G±, G 2 and H (see [16]). In order to simplify the 
notations, we assume that the small parameters of the long wave regime are equal (the general case can easily 
be obtained by modifying some constants) and set 

ei = C2 = H = e <C 1. 

The Boussinesq/Boussinesq system, in the one-dimensional case, can then be written using the convenient set 
of variables (771,772, wi,«2) = (Ci - C2,C2,d x ipi,d x ip2) as 

( d t m +d x (h 1 u 1 ) = 0, 
dtm +d x (h 2 u 2 ) = 0, 

dtui + d x {r]i + 772) + e(«i<9 x ui - -d 2 x d t u x - — dld t u 2 ) = 0, ( 4 ) 

1 + 3 6 
d t u 2 + d x (Ti]i + 772) + e(u 2 d x u 2 — — d 2 x d t u 2 - -d^dtUi) = 0, 

where hi = 1 + erji and h 2 = I + er/ 2 are the respective nondimensionalized depths of the upper and lower 
layer. Let us recall the following consistency result, that have been obtained in [16, Proposition 2.15] 

Proposition 2.1. Let U = (771,772, Ui, u 2 ) be a strong solution of system (1), bounded in L t ,oo ([0, T]; H s+to ) 
with s > 1 and to > 9/2, and such that (3) is satisfied. Then U satisfies (4) up to a residual R bounded by 



- * Mh 1 

\ firm 



\ R \ L oo H s < £ ^oIt— ,\U\ L i,oo H3+t0 

Remark 2.2. Here, and in the following, we denote by Co(Ai, A2, . . .) any positive constant, depending on the 
parameters Ai, A2, . . ., and whose dependence on Aj is assumed to be nondecreasing. Moreover, for < T < 00 
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and f(t,x) a function defined on [0,T] x R, we write / G L°°([0,T}; H s ) if / is uniformly (with respect to 
t G [0,T]) bounded in H s = H S (R) the L 2 -based Sobolev space. Finally, one has / e iv°°([0,T]; H s ) if 
/ G L cc {[0,T};H s ) and d t f G L°°([0,T]; H 8 ' 1 ). Their respective norm is denoted | • \ L ^ Hs and | • | L i,~ Hs . 

2.1. A new family of symmetric models 

In order to prove the well-posedness of hyperbolic systems as (4), one can use energy methods, that require 
symmetries of the system. Although our system is not entirely symmetrizable, we show in this section that it 
is equivalent (in the sense of consistency) at order 0(e 2 ) to a system of the form 

(So + e(Si(l7) - S 2 d 2 x )y t U + (s + e(Ei(tf) - ^ 2 dl))d x U = 0, (5) 

which satisfies the following crucial properties: 

Assumption 2.3. (1) The matrices So, £o, S 2 , £2 € Ai^M.) are symmetric. 

(2) Si(-) and £].(•) are linear mappings with values in .A/f^R), and for all U G M 4 , the matrices S\(U) and 
£i(J7) are symmetric. 

(3) So and S 2 are definite positive. 

Remark 2.4. We sometimes write the system (5) under the form 

p e (u,d x )d t u + Q e (u,d x )d x u = o, (6) 

with P t (U,d) = S Q + c(5i(l7) - 5 2 9 2 ) and Q t (U,d) = £ + e(Si(J7) - £ 2 9 2 ). 

First of all, let us point out that system (4) is only one among many other Boussincsq-type systems. In [5], 
Bona, Chen and Saut studied in the one-layer case a three-parameter family of Boussinesq systems, which are all 
approximations to the full Euler equations at the same order (in the sense of consistency) . The same structure 
applies also in the two-layer case, and we describe it quickly below. We then exhibit first order symmetrizers 
adapted to each of the Boussinesq-type systems, and leading to systems of the form (5). All of these models are 
equivalent in the sense that the original Boussinesq/Boussincsq system (4) is consistent at order C(e 2 ) with any 
of the systems presented in this section, as stated in Proposition 2.5. 

As a first step, one can use the following change of variables: 

v 2 = (1 - ea 2 9 2 )^ 1 w 2 , and v\ = (1 - e&i9 2 ) _1 (ui + eai<9 2 w 2 ), (7) 

with ai G M. and a2,&i > 0, and one recovers the three-parameter family of Boussinesq/Boussincsq systems 
introduced in [16]. 

Then, as it has been achieved in [5] , one can also use the classical BBM-trick [2] , and inherit new choices as 
parameters. This trick is based on the following calculations: since we have from (4) at first order 

d t m = -d x vi+0(e), dtm = -\d x v 2 + 0(e), ,. 

d t V! =-d x (rn+V2) + 0(e), d t v 2 = -d x { im + m ) + 0(e), W 

we can deduce the following, with the parameters Ai G [0,1] (i — 1 . . . 4), 

Sgvi = AiSgvi - (1 - Xi)d 2 x dtm + 0(e), 
dlv 2 = \ 2 3 x v 2 - 5(1 - X 2 )d 2 x d tV2 + 0(e), 

d 2 x d tVl = (1 - x 3 )d 2 x d tVl - x 3 d x ( Vl + 772) + 0(e), 

d 2 x d t v 2 = (1 - X 4 )d 2 x d t v 2 - X A dl{ lm + 772) + 0(e). 
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In the end, one obtains the following system, formally equivalent to (4) at order 0(e 2 ) 

d t U + A d x U + e(A(U)d x U - Aid%U - A 2 d 2 x d t U) = 0, 
denoting U = {r]i,r)2,Vi,V2) and 



A = 



(9) 



A 2 = 





/0 1 


M 






f ° 





Ai/?i 


-A 2 ai\ 








1 1 


1/6 



, Ai ■■ 





— A361 — A47ai — A361 — A 4 ai 







-A 2 — 
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\i i o 


0/ 


V-A47/32-A 3 i -A 4 /? 2 -A 3 i 





/ 


/(l-Ai)/3i 


(5(1 - A 2 )ai 
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/I* 


m o\ 






(1 - A 2 )a 2 




(1 - A 3 )6i (1 - A 4 )ai 


, ^(t^) = 






w 2 ?y 2 
fi 


V 
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1-A 3 )i (l-A 4 )/3 2 y 




\° 


u 2 / 



^ so that the system depends on 



with the parameters /3\ = 4 — 61, ct\ = ^ — 01, a 2 = gr^ — a 2 and j3 2 — a 2 -r j, 
the real parameters ai, a 2 , 61, and A^ (i = 1 . . .4), that can be chosen freely. 

In order to derive a system of the form (5), we exhibit a good symmctrizer of (9), that is to say 

S(U) = S + eS^U) - eS 2 dl, 

such that when we multiply (9) on the left by S(U), and withdrawing the 0(e 2 ) terms, we obtain a system (5) 
satisfying Assumption 2.3. 

The symmetrization of the one-layer shallow water system is well known, and consists in multiplying the 
velocity equation by the water depth. An adaptation of this to our two-layer model leads to 



So = 



/ 7 7 \ 

7 1 

7 

\ 1/5 J 



and 



Si(U) = 



/ 7U1 \ 

v 2 

7«i 7771 

\ v 2 m J 



(10) 



so that So and £0 = S0A0 are symmetric, and for all U £ R 4 , Si(U) and ^i(U) = Si(U)Aq + SqA(U) are 
symmetric. Moreover, So is definite positive for 5 > and 7 £ (0, 1): its eigenvalues are 



7, 1/5, and 1/2(1 + 7 ±v/(l- 7 ) 2 +4 7 2 ). 
Then, one can check that when we set 



/ a +(6- 1+7)6 

a + (56 

a+ (5- 1)6 + 7(61 A 3 - Ai/3i) 6 + 7A 4 ai 



V 



7 (||-A 2 ai) A4, V 2 " 2 



with a = 7((1 — A 2 )a 2 — (1 — Ai)/?i) and 6 = 7(1 — A 2 )cni, then 5 2 = SqA 2 + S 2 and E 2 = 5o-4i + S^o 
are symmetric. Then for any K £ R, we can substitute S2 + KSo for S 2 , and one has again S 2 and S 2 are 
symmetric. Moreover, since So is definite positive, one can choose K big enough for S 2 to be definite positive. 
Therefore, when we multiply (9) by S(U) = S + eSi(U) — e(S 2 + KSo)d 2 , and withdrawing the 0(e 2 ) terms, 
we obtain the perfectly symmetric system (5). 



Using the above calculations, it is now straightforward to obtain the following consistency result: 
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Proposition 2.5. Let U = (771, 772,^1,^2) be a strong solution of system (4) such that V = (771, 772, V\, V2), given 
by the change of variables (7), is uniformly bounded in L t '°°([0, T]; H s+5 ) with s > 1/2. Then V satisfies (9) 
and (5) up to a residuals bounded by € 2 C$ in the sense of L°°H S norm, with 

Co = Co(ai,a,2,bi,K,6+ t)\ v \ L ^,^ ^ 0T y H s+5y 

Proof. The first step in order to obtain (9) from (4) is to use the change of variables (7). Thus, when replacing 
in Tii by vi — eb\d 2 v\ — ea\d 2 V2 and U2 by V2 — ea2d 2 V2, we obtain straightforwardly that V satisfies (9) in the 
case 1 — Ai = 1 — A2 = A3 = A4 = 0, up to terms of the form 

e 2 dtd t f, e 2 d x (fd 2 x g) and e 3 d x ((d 2 f) 2 ), 

where / and g are components of V. Using the fact that V is bounded in L t :OO ([0,T); H s+5 ), and H S (M) is an 
algebra for s > 1/2, the remaining terms are clearly bounded by e 2 Co(ai, 02, b\, K, 5 + t> ) V , i,oo„ s+5 . 

In the same way, when we substitute the relations of the BBM trick (8) into the third-derivative terms of the 
equations, we obtain (9) up to extra terms bounded by e 2 Co(ai,a2, 61, j)|^Li.<»ij«+b- 

Finally, in order to obtain (5), we multiply (9) by So + eSi(V) — e(S 2 + KS )d 2 , and withdraw the terms 
e 2 (S 1 (V) - (S 2 + KS )d 2 x ) (A{V)d x V - AidlV - A 2 d t d 2 V). 
Each of these terms are clearly bounded by e 2 Co, since H S (R) is an algebra for s > 1/2, and 
\A(V)d x V\ HS < C \V\ LCOH3+3 , |5i(V)| ioo < C \V\ LOOH3 , and |^i^V| ff . + \A 2 d t d 2 x V\ Hs < C \v\ L ^ H3+1 

Therefore, the system (4) is consistent with the system (5) at the precision € 2 Cq. □ 

2.2. Well-posedness and convergence results 

The system (5) is a symmetric hyperbolic system, and can be studied using classical energy methods. We 
first prove in Proposition 2.6 that such a system is well-posed over times of order C(l/e). The proof uses in 
particular an a priori estimate of the solution, in an adapted norm: 

\U\ Hl+ , = \U\ H , + e\U\ H . +1 . 

Then, using a consistency result with energy estimates, we show that the solutions of our models converge 
towards bounded solutions of the full Euler system, assuming that such solutions exist (see Remark 1.2). 

Here and thereafter, we fix < 7 m i n < 7 < 7max < 1 and < S m - m < 5 < (5 max < +00. The limit cases 
(5 — > 0, 00 and 7 — > 0, 1) demand other scalings in the nondimensionalization than the ones presented in [16] 
(see Section A of [28] for example) and correspond to different regimes, such as the deep-water theory (from 
Benjamin [1] and Ono [35]), and lead to different models (see for example [7, 14] in the rigid lid configuration, 
and [26,36,39] in the free surface case). In all of these cases, the calculations of our justification break: 
the dependence of the constants Cq{ ,^_ ^ ,5 + ^) = Co(||So||, ||So" ||) m * ne following theorems prevents the 
parameters to approach the limits. 

In the same way, the following constants Co also depend on the set of parameters (ai, 02, b\,K,Xi) (i — 1 • • • 4). 
We decide to fix these to constants once for all, and do not write explicitly this dependency, in order to simplify 
the notations. 

We now state the well-posedness of our symmetric Boussinesq/Boussinesq model. 
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Proposition 2.6. Let U° € H s+1 , with s > 3/2. Then there exists a constant Co = Cb( rrz y ,fl + j) > 
suc/i t/iai /or e < eo = (CoK7° „ s+ i) _1 , there exists a time T > independent of e, and a unique solution 
U G C°([0, T/e);Hf +1 ) n ^([0, T/e);Hl) of the Cauchy problem (5) with U\ t=0 = U° . 
Moreover, one has the following estimate for t € [0,T/e]: 

\ U °\h° +1 

I^L-ao,*];^ 1 ) + \ dtU \L~(m;H°) ^ c °i-c \u°} , +1 et' (11) 

We postpone the somewhat technical proof to Appendix A. The key ingredients of the proof are quickly 
presented in the following remark. 

Remark 2.7. The condition s > 3/2 is necessary in order to obtain estimate (11), and thus the well-posedncss 
over times of order 0(l/e). One could obtain, using the exact same method as in the proof, the same result over 
times of order 0(1), with the sharper assumption s > 1/2 (which is the standard regularity for one-dimensional 
hyperbolic systems). 

The smallness condition (e|[/' |„ s+ i < Cq~ ) is also necessary, for the nonlinear terms in S\{U) to remain 
negligible when compared with Sq. Indeed, under this condition, the energy of the system (5), defined by 

E S (U) = 1/2(S A S U,A S U) + e/2(S 1 (U)A s U,A s U) + e/2{S 2 A s d x U,A s d x U), 

is uniformly equivalent to the • „,+i norm, that is to say there exists a > such that 

-Il7|l.+i < E S (U) < a\u\ 2 HB+1 . 
a < ' 

Moreover, the smallness condition on e \U°\ „ s+ i allows us to prove that the operator 

P e (U,d) = So + e(Si(U) - S 2 2 ) : H s+1 -»• H 3 ' 1 

is one-to-one and onto, and that P e (V ', 9) _1 Q e (V ', d) is uniformly bounded H* —> H 3 . This leads to 

\d t U\ H; = \P e {U,d)- 1 Q e (U,d)d x U\ H; < Co\U\ m+1 . 

Both of these properties are crucial in order to prove estimate (11). 

The existence and uniqueness of the solution of the Cauchy problem follow from the a priori estimate, using 
classical methods. 

Using the previous Proposition, and the consistency of the full Euler system (1) with our symmetric Boussi- 
nesq/Boussincsq model (5), one can now easily deduce the following convergence Proposition: 

Proposition 2.8. Let s > 3/2, e > and U = {Qi,C > 2,' i Pit^2) be a solution of the full Euler system (1) such 
thatV= {m,m,vi,v 2 ) e C°{[0;T/e); H° +1 ) n C^^T/e); H°) defined by 

V={Ci~ C2X2, (1 - eM*)- 1 ^^ + coi^wa), (1 - ea 2 d 2 x y 1 d x ^ 2 ) 

is uniformly bounded in L t '°°([0, T/e]; H s+5 ) and (3) is satisfied. We assume that e satisfies the smallness 
condition of Proposition 2.6, and denote by Vb the solution of the symmetric Boussinesq/Boussinesq system (5), 
with the same initial value Vbu =0 = VL =0 = V° and the same domain of existence. Then one has for all 
t€[0,T/e], 

\ V ~ VB \Loo([0,t];Hi+ 1 ) ^ £2fC< 0, 
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with Cq = Co(j^— , ,i_ y ,<5 + y, |V| 7 i,oo„ s+6 ,T). In particular, one has 



\V - V B \ LOO{[0tT/e] . H e+i ) < eC , 



with Cq independent of e. 



Proof. Let us first point out that the full Euler system (1) is consistent with (5) at the precision e 2 Co, that is 
to say that for any solution U of the full Euler system such that V G C°([0; T/e);H s+1 ) n ^([0; T/e); H s ) is 
uniformly bounded in L t '°°H s+5 ([0, T/e\), then V satisfies (5) up to a residual bounded by e 2 Co in the sense 
of L°°H S norm, with Co = Co(j^—,S + j, |V| ri,oo„, +6 ). Proposition 2.1 states that the full Euler system is 

consistent with (4) at the precision e 2 Co, and Proposition 2.5 achieves the result. 

Therefore, we know that V satisfies (5) up to e 2 /, with / e L°°([0, T/e];H s ), so that R s = A S V - A S V B , A s 
being the Fourier multiplier defined by A s u(() = (1 + \(\ 2 ) s / 2 u((), satisfies the system 

(So-eS 2 dl)d t R s +eA%S 1 (V)d t Ro+S 1 (Ro)d t V B ) + (j:o-e^2dl)d x R s +eA s (j: 1 (V)d x Ro+J:i(Ro)d x V B )^e 2 A s f, 

(12) 

With |/L«([0,T/ £ ];lf) ^ C o(^-^ + h\ V \ L l^H^- 

We can then carry on the calculations of Section A. 3, with the extra term e 2 A s /. We obtain that there exists 
Co( <^_ \ ,5 + g) such that as long as 

t\Ro\ H: +i < l/Co, (13) 

one has the estimate 

j f E(R s ) < eC (\V B \ H . + \V\ H .)\Ro\ 2 H . + e 2 (A s f,A s R s ), (14) 

with the energy E(R S ) defined by 

E(R S ) = l/2(S Rs,R s ) + e/2(S 1 {R )R s ,R s ) + e/2(S 2 d x R s ,d x R s ). 

Now, since So and S 2 are definite positive, the condition (13) implies in particular 

■^r\Ro\ H s+i < E(R S ) < Co | i?o | H s+i ■ (15) 

Go e e 

Thus, under this condition, and since V B is uniformly bounded in L°°([0,T/e]; H s ) with respect to e (from 
Proposition 2.6), one has 

d 
~dt 
and Gronwall-Bihari's Lemma leads to 



^-E(Rs) < eC E(R s ) + e 2 C \f\ Hs E(R s ) 1 / 2 , 



EiRs) 1 ' 2 < Coe\f\ Hs (e eCat -l). 

Finally, since -Roi t=0 = 0, and thanks to a continuity argument, there exists T(Co, /„,) > such that (13) 
holds for < t < T/e, and the estimate of the Proposition follows: 

\Ro\ H .+i < CoEiRs) 1 / 2 < Coe 2 \f\ Hs t. 

D 
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3. The KdV approximation 

In this section, we offer a rigorous justification of the so-called KdV approximation, from a class of symmetric 
systems that contains the symmetric Boussinesq/Boussinesq system (5). The KdV approximation consists in a 
decomposition of the flow into four parts, each of the components being lead by a Korteweg-de Vries equation. 
The construction of the KdV approximation is precisely explained in Section 3.1. Then, in Section 3.2, we 
obtain the convergence rate between the solutions of the coupled systems and the solutions defined by the KdV 
approximation. As a consequence, when we combine this result with the convergence Proposition 2.8, it follows 
immediately that any strong solution of the full Euler system existing over times C(l/e) and bounded in a 
sufficiently high Sobolev norm, is well approximated by the KdV approximation. More precisely, we state the 
following: 

Theorem 3.1. Let s > 3/2 and U — (Ci, C2, ipi, 1P2) be a solution of the full Euler system (1) such that 
V = fa, r)2,Vi,V2) e C°{[0;T/e);H s e +1 ) nC'dO;^);^) defined by 

V = ( Ci - C2 , C2 , (1 - ehSft-Hd^i + e ai d 2 x v 2 ) , (1 - ea 2 d 2 x )- 1 d x ^ 2 ) 

is uniformly bounded in L t '°°([0,T/e]; H s+5 ) and (3) is satisfied. Then there exists 1 (ei, . . . , e^) a basis ofM 4 , 
and coefficients Ci, Aj, Hi (i = 1 • • • A) such that, denoting by m the solution of the KdV equation 

d t Ui + Cid x Ui + eXiUidxUt + e^d^Ui = (16) 

with u»l =0 = 6j • 5oV[ t=0 (So defined in (10),), one has for all t € [0, T/e], 

4 

\v -Y J u ^\ L ~([o.t];m +1 ) ^ e ^~ tc °> 

with C = Co{j^-, 7 (i_ 7 ) ^ + h\ V \ L ]-°°H°+^- 

Moreover, ifV satisfies (1 +x 2 )V| t=0 € H s+ , then one has the better estimate 

4 
\V -2_^u l e l \ LOO([0T/<iyH;+1) <eCo, 

i=l 

with C - C (^-, 7T i L 77 ,'5 + h \V\ L] ,oo H3+5 , |(1 + x 2 )V u J Hs+i ). 

Proof. The proof proceeds from different results of the paper; the completion is as follows. In Proposition 3.12, 
we prove the convergence between the solutions of systems of the form (17) (thus containing symmetric Boussi- 
nesq/Boussinesq systems (5)) and the approximate solution U app = Yli=i u i e i + AJ\, defined in Definition 3.6. 
The residual U\ is then estimated in Proposition 3.9. Finally, since we have from Proposition 2.8 the convergence 
between the solutions of the full Euler system (1) and the solutions of the symmetric Boussinesq/Boussinesq 
system (5) with a better rate, the Theorem follows; see also Remark 3.13 below. □ 

Remark 3.2. The difference on the convergence rate for different sets of initial values is not simply a technical 
issue. Indeed, one can see in Figure 7, page 30, that if the condition of sufficient decreasing in space is not 
satisfied, the convergence will be worse than 0(e). Requiring the initial data (and thus the solutions of the 
KdV equations) to lie in weighted Sobolev spaces guarantees that the nonlinear interaction between the four 
traveling waves can be neglected. As a matter of fact, this condition on the sufficient decay in space of the 
initial data appears also naturally for the KdV approximation of the one-layer problem, as we see in [6,38]. 



one has explicit expressions for the basis (ei, . . . , 64) and the coefficients c;, A;, ^t;, that we display in Remark 3.5 page 13. 
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3.1. Formal derivation 

The class of system that we now study is the following: 

(s + e(Si(E0 - S 2 d 2 x ))d t U + (e + e(E!(t/) - E 2 ^))d x t/ = 0, (17) 

with the following hypothesis: 

Assumption 3.3. (1) The matrices So, So, S 2 , E 2 € M^M.) arc symmetric. 

(2) Si(-) and Ei(-) are linear mappings with values in Al^R), and for all U € R 4 , S±(U) and Ei(/7) are 
symmetric. 

(3) So is dchnite positive, and Sq So has four different non zero eigenvalues c, (i = 1 . . .4). 

Remark 3.4. The symmetric Boussinesq/Boussinesq systems (5) derived in Section 2.1 immediately satisfy 
Assumption 3.3, with 



Ci_± V 26 

Following the classical WKB method, we look for an approximate solution of the Cauchy problem (17) with 
initial data U under the form 

U app (t,x) = U (et,t,x) +eUi(et,t,x), 
with the profiles Uo(r,t,x) and dJ\{T,t,x) satisfying Uo\ t=T=0 = U° ct Ui\ t=T=0 = 0. 

We plug the Ansatz into (17), and obtain 

(S d t + Z d x )U + eSod T Uo + e(Si(l7o)$E>b + Zi(U )d x U - S 2 d 2 x d t U - E 2 ^t/ ) 

+e(S d t + E a x )C7! + e 2 R = 0. (18) 

We now deduce the equations satisfied by Uo(r,t,x) and Ui(r,t,x), solving (18) at each order. 

At order 0(1): We solve 

(S d t + 2od x )U = 0. (19) 

We assumed that So is symmetric definite positive, and hence induce a scalar product: 

(u, v) = u ■ Sol) = u SqV. 

Since Sq 'So is real and symmetric for the scalar product (•, •), it is diagonalizable in an orthonormal basis. We 
denote by e, (i = 1 . . .4) the basis vectors, which are the unitary eigenvectors of S^Eo- By definition, they 
satisfy for 1 < i, j < 4, 

(ei,So~ E ej) = ej • E ej = c l 5 l , J , and (e^e,-) = e t ■ S ej = Sij, 

4 
with Si.j the classical Kroneckcr delta symbol. Therefore, when we define ui = ei-SoUo (and hence Up = Y] u^e,), 

multiplying (19) on the left by e*, we obtain 

(d t + Cid x )ui = 
for all i = 1, . . . , 4. Finally, since Uj satisfies a transport equation, we use the notation 

Ui(r,t,x) = m(T,x-Cit) = Ui(r,Xi), (20) 

with initial data Uj(0, Xi) = e« • SoU°(xi). 
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At order 0(e): We solve 

S d T Uo + £i(t/o)dz^o + Si(U )8tU - Z 2 d 3 x U - S 2 d 2 x d t U + (S d t + S 5 x )f/i = 0, 

that we can split 2 in 

d T m + X l u l d Xz u i + Hidl.m — 0, 

with Xi = e, • (Si(ej) — Cj/Si(ei))e, and /Uj = ej • ( — £2 + CiS2)ei] and in the other hand, 
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^ a ijk u k (T, x - c k t)d x Uj(T, x - Cjt) 
(j,fe)^(*>») 



with a,' 



(d t + c» ck )(ej,i7i) + 
= e, • (Si(efe) - c J -5 , i(e fc ))e j and /3,j = ej ■ (£ 2 - CjS 2 )e 



^^ y a|uj(r,a;-Cjt), 



(21) 
(22) 

(23) 



It is clear that «j satisfies (20) and (22), if and only if m(et,t,x) satisfies the Korteweg-de Vries equation of 
Theorem 3.1: 

d t Ui + Cid x Ui + e\iUid x Ui + enid^m = 0. (24) 

Remark 3.5. In the specific case of symmetric Boussinesq/Boussinesq systems (5), one obtains the following 
values for the coefficients: 



(ci)ie{i...4> = (c+,-c+,c-,-c-), 
(Ai)ig{i...4} = (A+, A+, A_, A_), 

(Mi)»e{i...4} = {n+,-H+,H-,-H-), 



with c ± = 



with A 



± 



with /i± 



1 + ^^(1-5)2 


+ 47(5 




2(5 

3 (2 - 5)01 + (5 - 


l-(l-7) 




2 e±( c+ - 

c± (! + £ + £)( 


4-fc?)- 


±r 2 

<5 C ± 



c ± / (5+1 



with 0± 



2<5(c 2 - c 2 



1 1 . The unitary eigenvectors of S £q are given by 



ei 



1 

e7 



i 



V«5c + -(5/ 



, e 2 



e. 



/\-i \ 



\Scl-Sj 



e.j 



e_ 



, e 4 



e. 



\5-Sc 2 _J 



( ±A 



\S-5c 2 _J 



Let us first remark that these coefficients only depend on 7 and 5, so that neither the change of variables (7) 
nor the BBM-trick (8) affect the coefficients of the KdV approximation. 

It is also worth pointing out that the dispersion coefficients /j,± cannot be zeros for 7 G (0, 1) and 5 > 0. 
Therefore, the KdV approximation cannot degenerate into Burgers-type equations. 

The coefficients of the KdV equation as an approximation to describe solitary internal waves in the free 
surface case have already been introduced in [22], and then in [30,32,34], and correspond to the ones we present 
here. In a slightly different regime, Craig, Guyenne and Kalisch [14] obtained a model that consists in four 
independent propagating waves, two of them satisfying the KdV equations of our slow mode waves (the other 
two being solutions of the first order transport equation d t u ± c + d x u = 0). We state here that in the long 
wave regime, any bounded solution of the full Euler system can be decomposed into four propagating waves, 
each of them being well approximated by independent solutions of four KdV equations. This simultaneous 
decomposition of the flow is new, and allows to compare with other models, such as Boussinesq-type models. 



This splitting is in fact necessary. Indeed, the multiscale WKB expansion can be justified for times of order 0(1/ 1) only if the 
growth of the corrector term Ui is sublincar. As wc sec in Proposition 3.9, thanks to Lemma 3.8, the particular form of (23) allows 
to obtain a square-root growth, and even better if U° is sufficiently decreasing in space. 
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3.2. Rigorous demonstration 

The strategy is the following. Using the previous calculations, we define the approximate solution: 

Definition 3.6. Let U° <G H S+2 (M), with s > 1/2. We call approximate solution of the system (17) any function 
U app {t,x) e C°([0,T/e);H s ) 4 such that 



U app (t,x) = ^u l (t,x)e i + eUi(et,t, 



(25) 



where (u,)j=i...4 is the solution of the four uncoupled KdV equations (24), with Ui\ t=0 — e^ • SqII , and the 
correcting term U\ is a solution of (23) with U\\ T=t=a = 0. 

We first prove that such solutions exist in the strong sense for sufficiently smooth initial data, over times of 
order C(l/e). Then we use estimates on Uq and Ui, to obtain a consistency result. This result allows us to show 
that Ua,p P (t, x) indeed approximates the solution of (17) with the same initial data, at least at order 0(e 3 ^ 2 t). 

Proposition 3.7. Let U° € H S+2 (R), with s > 1/2. Then, one has: 

(1) For all i = 1 . . . A, there exists a unique strong solution of the Cauchy problem (22), with initial data 
u i\ t=0 — e i ' SqII . Moreover, there exists T > 0, such that one has the estimate 



E 



" ,l I»([0,T];iI^)- C '»\ 7 ( 1 _ 7 )'- ' g 



1 



1 



\u°\ 



ffs + 2' 



(2) There exists a function U± € L°°([0,T] x M;_ff s ), strong solution of (23), with Uu = 0. 

Proof. (1) The existence and uniqueness of the solutions of (22) is classical: see [24] for the local well- 

posedncss of the KdV equation, and [12] for the global well-posedness. One obtains the estimate as 
usual: as we multiply the equation (22) by A 2k Ui (with 3/2 < k < s + 2) and integrate with respect to 
the space variable, one obtains 



II ( (K k Ul ) 2 dx 
2 dt J R V l> 



\i I K k {uid x Ui)h k Uidx 



Thanks to the Kato-Ponce Lemma, we manage to estimate the right-hand side as follows: 



A k (u i d x u l )A k u l dx 



< 



fc„, ^2 



d x Ui(k k Ui) 



[A k ,u,i\d x u t (A k Ui)dx 



< Co u 



»|JT*' 



and we conclude by applying Gronwall-Bihari's lemma, that reads 

Ui I < r Ku=o|gfc 

\Ui\ Hk S Oo- — —, J — , 

1 - C a t\u lU=0 \ Hk 

so that the estimate of the proof follows for T sufficiently small, since 



\ui\ t = \ H3+2 = |e,-S* C/ \ Hs+2 <C y^_ yS+fJ \ u Ihs+2- 
Let us recall the notation (20): we have thus proved that m(T,t,x) <G L°°([0, T] x R;7J S+2 ), with 

4 



Z^ \ Ui \L°°(lO,T]xR;H°+ 2 ) - C ° I 7 (- 1 _ 7 ~ ) '' 5 + A ) \ U I// -" 



1 



TITLE WILL BE SET BY THE PUBLISHER 15 



(2) We can then exhibit XJ\\ let us write (23) under the form 

(d t + Cid x )(ei, Ui) = ^2 fijk{r, t,x)+J^ j d x g^{r, x - Cjt). 

{j,k)jt(i t i) j^i 

From the above estimate on m, one has fijk G L°°([0,T] x K;iJ s+1 ), and gij G L°°([0,T]; H s ). Hence, 
for s > 1/2, one can set 

/"* 1 

(e i ,U 1 )(T,t,x)= ^2 / fjjkJT, s,x + Cj(s - t))ds + y^ — - — (gij(T,x-Cjt)- gij(r,x-Cit)), 

and U\ satisfies the hypotheses of the Proposition: Ui(r,t, x) G L°°([0,T] x R;H S ) and J/n =t=0 = 0. 

□ 

We now prove that Ui, which is the corrector term defined by (23), and that contains all the coupling effects 
between the different components, obeys to a sublinear secular growth. The key point is given by the following 
Lemma, that proceeds from Propositions 3.2 and 3.5 of [29]: 

Lemma 3.8. Let u be the solution of 



{d t + cd x )u = g{v u v 2 ) mth Vze{12} ( (d t + c t d x )v l= 



with c\ ^ C2, v\, v\ G H S (M), s > 1/2, and g is a bilinear mapping defined on M 2 and with values in M. Then 
one has the following estimates: 

(1) //c = ci , then lim j t \u(t,-)\ He{R) ^0. 

(2) Ifc^ Cl ^c 2 , then±\u(t,-)\ Hs{R) =0(l). 

Moreover, if there exists a > 1/2 such that vf(l + x 2 ) a , and V^O- + x 2 ) a G H S (W), then one has the better 
estimate 

\ u \l^H'(r) - Co\vi(l + x ) a \ H3{R) \v 2 (l + x )"| ff3(R) , 

with Cq = Co(c, ci, C2). 

We are now able to give the following crucial estimate on U\. 

Proposition 3.9. Let s > 1/2 and U° G H s + 2 . Then with U x G L°°([0,T] x R;H S ) a strong solution of (23) 
with U\\ „ „ = 0, one has the estimate: 

J -|-r = t = ' 

\ Ul \L°°([0,T]xlO,t];H>) - CaV ^> 

With Co = C (^y,5+ \, \U°\ Hs+2 ). 

Moreover, ifU satisfies U (1 + x 2 ) G H S+1 (R), then one has the uniform estimate 

W Al<*>{[0,T]xWL;H>) — ^0\U (l + X )\ H . +1 , 

with Co = Co{^ T) ,5+ |, \U°\ H3+2 ). 
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Proof. Let us decompose U\ as a sum of functions as in Proposition 3.7: 



/"' 1 

(ei,Ui)(T,t,x) = J^ / /jjfc(T, s,x + Cj(g - t))rfs + y^ — - — (gij(r,x - cjt) - (gij(r,x - Cjt)) 

(j,fe)#(i,i) - 70 j^. Ci Cj 

with the functions /^fc and g^- coming from (23) written in a simplified form: 

(<9 t + c l d x )(e l , Ui) = ^ fijkir, t,x) + ^ d x9ij(T, x - c 3 t). 

(j.k)^(iS) &i 

From Proposition 3.7, we know that the following bounds hold: 

Vt€[0,T], 2^ |/y fc L=»([0,T]xM;H*+ 1 ) + Z^ \9n\L°°{[Q,T\\H>) - C °\ U \h°+ 2 > 



with Co = ^0(^^,5 + ^). 



Therefore, one has 



1 „ 1. 



Moreover, again for j ^ i, we remark that fijj can be written as 

so that C/ 4JJ ' has the same form as V 4 -?, and can be treated in the same way. And since /^j e L°°([0, T]xR; H s+1 ), 
U 1 ™ is uniformly bounded in i/ s . The last terms that have to be bounded are U lJ , for (j, fc) 7^ (i,i) with 
j 7^ fc, [7* J . One can easily check that XJ % i k satisfies the hypothesis of Lemma 3.8, with fijk — g(uj,d x Uk), for 
any r € [0, T]. We then immediately deduce 

\ U AL^([0,TMO,t];Hs ) <^ C o( lil _ i yS+y\U \ H3+2 ). 

As for the second estimate of the proposition, let us first remark that the estimates of V^ and IP™ are 
timc-independant, and in agreement with the improved estimate. Therefore, the only remaining terms we have 
to control are U^ with j 7^ k. Of course, we will use the second case of Lemma 3.8, but we have to check first 
that for every r € [0,T], the initial data Uj(t, 0,x) and d x Uk(T 7 0,x) are localized in space, that is 



Vre[0,T], \(1 + x 2 )u 3 (t,0,x)\ Hs + |(1 + x 2 )d x u k (T,0,x)\ Hs < 



do. 



This property is true at r = (by hypothesis of the proposition), and is propagated to t > 0, using the fact 
that Ui(r, Xi) satisfies the KdV equation (22). This propagation of the localization in space has been proved by 
Schneider and Wayne in [38, Lemma 6.4]. We do not recall the proof here, and use directly the statement: 

Lemma 3.10. If (1 + x 2 )UP_ a £ H s+1 , then there exists C\,C\ > such that 

|(l + x 2 K(r,0,x)| HS+1 < C x I (1 + aJ 2 )«j k _ t _ 1^. +1 < Ci|(l + x 2 )£/°| t=0 | HS+1 . 
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This Lemma, together with the second estimate of Lemma 3.8, allows to control W^, uniformly in time. 
Every term of the decomposition of U\ has been controlled, and one has the following estimate: 

|tfiL~ ([0 ,:n xR;ff .)< C °(^i^y S+ - 5 >\ U0 \Hs + i)\(l + x 2 W \ 2 H s +l - (27) 

This concludes the proof. □ 

The next step consists in proving the consistency of our approximation with the symmetric system (17). 

Proposition 3.11. IfU € H s+5 with s > 1/2, then U app (t,x) defined in Definition 3.6 satisfies the symmetric 
system (17) up to a residual of order 0(e 3 / 2 ) in L°°([0, T/e]; H s ). 

Moreover, Uq satisfies U°(l + x 2 ) G H S+4 (M.), then the residual is uniformly bounded in L°°([0,T/e\; H s ) by 

e 2 C Q , with Co = Co (^),S + J, |(1 + x 2 )V°\ H . +t ,T). 

Proof. Plugging f/ app (i, x) into (18), we see from the calculations of Section 3.1 that the only remaining term 
we have to control is e 2 R(et,t,x), with 

R = d T Ui + Si^a^i + £i(I/i)0*E/b + S^Uo)^! + Si(J7i)ftl7o - £2^1 - S 2 d 2 d t U 1 

+eE 1 ((7 1 )^(7 1 +e5 1 (C/ 1 )a x ;7 1 , 

where U (et,t,x) = X)»=i M»(*,a;)ei. 

Each term of the right hand side is suitably bounded in the Sobolev iJ s -norm, as we show in the following. 
Indeed, from Proposition 3.9, one has immediately 

\22%U 1 (et,t,-)\ H .<Co\U 1 (et,t,-)\ H . +a < VtC (-^-^- y S +^\U \ hs+ ,,tY 

Then, from (23), we deduce 

(e l ,d t U 1 ) = -Ci{ei,d x Ui) + fi 
with fi e i°°([0,T] x M;i7 s+2 ), and |/ 4 | H3+2 < C Q \U°\ He+5 , so that one has identically 



I^M&EMeM,-)!*. <Co|9 t £/i(ef,t,-)| g . +2 <VfCoC t^ ■ ■" "~ S r " 



1 

7(1-7)'" ' £' r lff3+5 

One obtains in the same way the desired estimates for Y>i(Uo)d x Ui, J^i(Ui)d x Uo, Si(Uo)dtU\, S\(Ui)dtUo, 
Zi(Ui)d x Ui and Si(Ui)d x Ui. 

Finally, in order to estimate d T Ui, we differentiate (23) with respect to r. Since m satisfies (24), one has 
d T Ui € L°°([0,T); H s+2 ). We are on the frame of the Lemma 3.8, so that we can obtain as in Proposition 3.9 
that drUx e i°°([0,T] x [0,t];H s ), with 

\d T U x {cb,t,-)\ H . < C Vt\d t U°\ HS+2 < C Vt\U°\ HS+5 . 

Hence, R(et,t,-) e L°°([0,T/e]; H s ), and 

\R\h' — Co\/T/e\U \ HS+5 , 

which concludes the first part of the proof. 

The second part follows in the exact same way, using the second estimate of Proposition 3.9. □ 

Finally, thanks to the consistency result and the estimate on U\ , we are able to set the following convergence 
Proposition: 
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Proposition 3.12. Let U° e H s+5 , s > 1/2, U B e L co ([0,T/e];H s + 5 ) be a family of solutions of (17) with 
Ub\ 1=0 = U° and U app be defined by Definition 3.6, with the same initial value. Then one has 

\U app - u B\ Laa([oA . H; +i ) < C e 3/2 t, 

with Co = Co(^ y ,<5+ i, \U°\ H3+5 ,T). 

Moreover, U° satisfies U°(l + x 2 ) G _ff s+4 (M), then one has the better estimate 

\U ap p - C/ B| L oo([ ,t];ff|+ 1 ) - C 0C t, 

with Co = C (^ ) ,6+ \, \U°(l+x 2 )\ Hs+i ,T). 

Proof. Let us set R c = U app — Us- Thanks to Proposition 3.11, we know that 

(So - eS 2 d 2 x + eS^U^OtR" + (S - eS 2 9 2 + e£i(f/ ap p))d x .R e = e^ 2 f + eA + eB, (28) 

with A = d t Si(U avp )R e - SxiR^dtUs, B = 9 x Si(J7 app )i2 € - ^(R^OJIb and a function / € L°°H S . 

Then, we can follow the same path as for the proof of Proposition 2.6, in Section A. 3 (see also the proof of 
Proposition 2.8). We define the energy as 

E S (R<) = ±(S A s R e ,A s R e ) + e -{S 2 d x k s R\d x k s R<) + |(^(C/ a pp)A s i? e , A s i? £ ), 
and the exact same calculations lead to the following inequality: 

^-E s {tf) < CoeE s {tf) + Coe*l 2 {E s {R*)) 1 ' 2 , 
at 

with Co = C (^ T) ,6+ i, \U°\ HS+5 ). 

From Gronwall-Bihari's theorem, we get E s (R e ) < Coe 1 l 2 {e Catt - 1), and finally for et < T , 

\V app -U B \ H . + i <CoE s {W)<Co{ \ 8 + - 5 ,\U°\ Hs+& ,T)e z ' 2 t. 

The second part of the proof follows in the same way, using the consistency at order 0(e 2 ) of Proposition 3.11. 

□ 

Remark 3.13. From this Proposition, one can immediately deduce the convergence rate between the solution 
of the symmetric system (17), and the KdV approximation (24). Indeed, since we know from Proposition 3.9 
the growth of the correcting term U±, and since we restrict ourselves to times < t < T/e, the convergence rate 
is of order O(evt) in general, and of order 0(e) if the initial data is sufficiently decreasing in space. 

In the same way, we obtain the convergence rate between bounded solutions of the full Euler system (1), and 
the KdV approximation (24), using Proposition 2.8. This result is stated rigorously in Theorem 3.1. One sees 
that even in the case where the initial data is rapidly decreasing in space, the convergence rate of the symmetric 
Boussincsq/Boussincsq model (5) (namely 0(e 2 t)) is better than the one of the KdV approximation (namely 
0(e)). This is due to the interaction between the traveling waves of different wave modes, that is captured by 
the Boussincsq/Boussincsq system, and not by the uncoupled KdV approximation, and which is of order 0(e) 
for times of order 0(1). The decreasing in space of the initial data allows this error to remain of order 0(e) for 
times of order 0(l/e). 

Numerical simulations for both the Boussincsq/Boussinesq models and the KdV approximation are computed 
in Section 4.2. In particular, the relationship between the convergence rate and the decreasing in space of the 
initial data is discussed and enhanced in Figure 7. 
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3.3. The models under the rigid lid assumption 

In this section, we formally recover models existing in the literature in the rigid lid configuration. Start- 
ing from our Boussinesq/Boussincsq model (4), we recover the three- parameter family of rigid lid Boussi- 
nesq/Boussinesq systems presented in [7] . One can then apply the method presented in the previous section, in 
order to obtain the KdV approximation in this case. 

The rigid lid models use the variables (C,v), where £ is the interface deviation (— r)i = r\ 2 = C)i an d v is the 
shear velocity defined by 

v = (d x (/)2 - r yd x (j)i)\ z=eC 
Using the calculations in [16], one has 



v = u 2 - 7Wl - e [^diu, + ^— + ^ j d z x u 2 ) + 0(e 2 ). 

Then, adding the first two equations of (4) leads to d x {h\U\ + h 2 u 2 ) = 0, so that one has u\ + \u 2 = 0(e), and 

v = 5 -^ 1 u 2 + O(e) = -(5 + 7)«! + O(e). 


Therefore, using a straightforward combination the equations of (4), one checks that the system becomes 



dtC + s^B x v + e^d x ((v) + e^$-dlv = G(e 2 ), 
d t v + (1 - 7 )0 X C + efefrvdxv = 0(e 2 ). 



(29) 



Finally, using BBM- tricks as in (8), and the change of variable vp = (1 — e/3<9 2 ) 1 w (with (3 > 0), one obtains 
eventually the three-parameter family of rigid lid Boussincsq/Boussinesq systems presented in [7] : 



(1 - ebd 2 x )d t C + j^d x vp + e ^d x ((v ) + mdlvp = G( £ 2 ), 

(1 - eddl)d t vp + (1 - 7 )5xC + t-$^iVf}d x Vfi + ec(l - 7 )^C = 0(e 2 ), 



(30) 



with a, b, c and d set (with 6\ > 0, 62 < 1, /3 > 0) as 

From this system, one can easily follow the path of Section 3.1, and deduce the KdV approximation related 
to system (30). One would obtain a similar result as in Theorem 3.1. Eventually, the KdV approximation 
consists in decomposing the approximate solution (CKdVj^Kdv) as 

(CKdv, w K dv) = u+e + + u_e_, 

with u+ and U- two solutions of independent Korteweg-de Vries equations, namely 

dfU± ± cd x u± + e\±u±d x u± ± e/j,d x u± = 0, (31) 

with the following vectors and coefficients: 

1 f±-^=\ fl^ . 1 3c <S 2 - 7 c(l + 7 <5) 



y/2 \VjTSj ' ' \j + 6' ^2(1-7) 2 7 + ^ 6 <5( 7 + 5) ' 



20 TITLE WILL BE SET BY THE PUBLISHER 

In that way, when looking at the decomposition of deformation of the interface, the KdV approximation 
leads to two counter-propagating waves, that is to say that one can write r\ = 77+ + rj—, with r\± solution of 



dtV± ± yTT5 dxfl± + e YTT5 V±dxT1± ± e 66fr+S)'' h: = 

We recover the classical KdV equations in the rigid lid configuration (see for example [14,15,26,33]). Following 
Section 3.2, one would obtain in the same way a rigorous justification for the KdV approximation, under the 
rigid lid assumption. 

One sees that whereas the KdV approximation in the rigid lid case leads to a decomposition into two waves, 
the free surface configuration predicts the decomposition into four waves, each of them solution of a KdV 
equation with different velocities. This striking fact leads to think that the rigid lid assumption may induce a 
significant alteration of the behavior of the solutions, and thus cannot be considered as a harmless statement in 
all configurations. This remark has already been addressed in [25,34], but to our knowledge, in the absence of 
the exhaustive decomposition given in Theorem 3.1, the analysis of the difference between the two configurations 
has never been extensively discussed. 

The following Section is devoted to a detailed study of the differences between the rigid lid and free surface 
configurations, depending on the values of the density ratio 7 and the depth ratio 5. 

3.4. Discussion 

Let us recall here that the KdV approximation in the free surface configuration consists in decomposing 
U = {rii,ri2,Ui, U2) as U ~ Yli=i u i e ii with m satisfying the KdV equation 

d t Ui + c % d x Ui + eXuidxiit + e/U,d 3 Uj = 0. 

The coefficients c;,Aj and /jj, as well as the vectors e^ are given in Remark 3.5, page 13. This leads to the 
following decomposition for the respective deformations of the interface and the surface: 

4 
(1) For the interface: r/2 = / Ui^i.2 = / Vj,k, and r]± t ± satisfies the KdV equation 

4=1 (j,fc)=(±,±) 

d t ri±,k ± c fe 9 x 77±,fe + e\ l k r]±,kd x r] ±tk ± e^fc^»7±,fc = 0, 

where ejj denotes the j th component of the vector ej and with the following coefficients 3 : 

2 _ l + 5±y/(l-6)*+4r,6 _c± (l + ^ + ^)(4-j+j £ )-H 

' " 26 ' M± 6 d -2^ 



X\ 



3 c± (2-<5)4 + <5-i-(l- 7 ) 



■ + 

4 



\(cl-d)(l-4 



(2) For the surface: (1 = r}\ +772 = \. u i( e i,i + e «.2) = /^ Cj.k, and C±,± satisfies the KdV equation 

»=i (j,fe)=(±,±) 

d t (±,k ± c±d x (±, k + e\ s k C±,kd x C±,k ± ^kdlC±, k = 0, 

with the same values as previously for c± and //±, and 

_ 3c ± (2 - ,5)4 + S - I - (1 - 7) 
A+ — — — 



^± 



(4-c 2 _)c^ 



These are the coefBcients that are displayed in [22, 30, 34] . 
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Figure 2. The different velocities, in the rigid lid configuration (c = ./-4) and free surface 



case 



(c± = yi7^7(pw ); for 7 G (0) 1} and (a) s = 1/2) (b) 5 = 2 



Therefore, one sees that for both the surface and the interface elevations, the KdV approximation predicts the 
evolution of four different waves, two of them corresponding to the velocities ±c+ (we call them fast mode waves), 
and the other two corresponding to the velocities ±c_ (we call them slow mode waves), with c+ > c- > 0. 
The fact that such different modes exist is characteristic of the free surface configuration, as only two counter- 
propagating waves appear in the rigid lid case. Moreover, as we see in Figure 2, the two velocities corresponding 

to the free surface configuration can be very different from the velocity in the rigid lid case (namely c = WtEtJ), 
depending on the values of 5 and 7. In these cases, one expect the solutions in the two configurations to behave 
very differently. 

The aim of this section is to study more in depth the behavior of the KdV approximation in the two different 
configurations, with respect to the parameters 7 and S. The first part is devoted to the study of solitary waves, 
as known solutions of the KdV equations. Then, we study the case where the initial data satisfy the rigid 
lid hypothesis, and explore the evolution the surface in that case. The results we obtain are summarized in 
Section 3.4.3; we let the reader refer to Figures 8-15 for a numerical illustration of our statements. 

Remark 3.14. In the following study, and especially in Table 1, we allow ourselves to look at the behavior 
of the system in the limit cases of the parameters (5 — > 0, 00 and 7 — > 0, 1), despite the fact that the rigorous 
justification of the KdV approximation, as well as the Boussinesq/Boussinesq models, break in these limits. 
However, the KdV approximation when 7 — >■ 1 has been widely used in the literature (see for example [19,21, 
37,39]), and such limits offer striking illustrations of our discussion. 

3.4.1. Solitary waves 

It is well known that the solitary wave solutions of the generic KdV equation 



dtu + cd x u + eXud x u + e/id x u = 



can be expressed as follows: 



u(t, x) 



M 



cosh(fc(x — xq — c't))" 2 



(32) 



with c' = c + e^jr-, k = \I j^-, and M and xq arbitrary. 
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We discuss in the following the polarity, magnitude and thickness of such waves, as the parameters 7 and 8 
specify the coefficients of the KdV approximation. 



Polarity. It is obvious that for k = \/j^- to be real valued, the sign of the ratio X/fj, determines the sign of the 

acceptable values of M . Therefore, we are able to predict, depending on the parameters 7 and 6, the polarity 
of the solitary waves predicted by the KdV approximation (elevation or depression). We give here the result, 
first for the free surface case, and then in the rigid lid configuration. 

First, one can check that for every value of 8 > and 7 £ (0, 1), the three coefficients A+, A+ and fi + > 
are positive. Hence, the fast mode solitary waves will always be of elevation/elevation type (both surface wave 
and interface wave are convex upward). The qualitative nature of the fast mode is thus similar to that of the 
one-layer water-wave problem, which is always of elevation type. In particular, when we set 7^0 + and <5 >■ 1 , 



one recovers the classical KdV equation for a single layer at the interface (c + 



1, X* 



3/2, M+ -> 1/6), 



and when we set 5 — >• 00, one recovers the classical KdV equation for a single layer at the surface (c+ — > 1, 
X a + -> 3/2, n + -> 1/6). 

The behavior of the slow mode is more peculiar, as the nonlinear coefficients X s ! 1 can have both signs, 
depending on the size of the thickness ratio 8. Indeed, as it has been pointed out in [22,41] and then in [34,36], 
for every 7 <E (0,1), there exists a critical ratio <5 C (7) such that if 8 > 8 C (^), then X l _ > and Af_ < (and 
conversely if 8 < 8 C ). Since one has /i_ > for every value of 8 > and 7 £ (0, 1), we know that the slow mode 
solitary wave will be of elevation/depression type (surface wave convex upward, and interface wave concave) if 
8 < 8 C , and of depression/elevation type if 8 > 8 C . There is no solitary waves in the case 8 = S c . 

More precisely, the critical ratio is the unique real solution of the equation 

X 3 + ( 7 2 + 37 - 3)AT 2 + (3 - 4 7 )V -1 = 0, 



and takes values in 8 C £ (1, 5/4] for 7 £ (0, 1) (see Figure 3). 

In that way, the behavior of the slow mode waves 
resembles that of the interface waves with a rigid lid. 
Indeed, such a critical ratio appears straightforwardly 
in the rigid lid configuration: if 8 < 8' c = y/j, then the 
interface wave is of depression type, and if 8 > 8' c , then 
the interface wave is of elevation type. This critical 
ratio is well known in the literature, and has led to 
many extended models, where a cubic nonlinear term 
becomes the major source of nonlinearity for 8 ~ 8 C 
(see for example [15,18,20,21,26]). 

It is interesting to see that even if the respective 
polarity of the interface slow mode waves in the free 
surface configuration, and the interface waves in the 
rigid lid configuration follow qualitatively the same 
behavior, the value of the critical ratio is notably dif- 
ferent. Indeed S c is always located above 1 in the 
free surface case, on the contrary to the rigid lid case. 
Moreover, as we can see in Figure 3, they have con- 
siderably different values, except when 7 ~ 1. In the 
area between the two curves, the polarities of the in- 
terface waves in the free surface and in the rigid lid 
configurations are reversed; the two models therefore 
lead to considerably different results. 



1.4 



1.2 



1 • 



0.8 



0.6 



0.4 



0.2 



r—--___ 


• 
/ 




- / 
/ 

1 


free surface configuration" 

■ - rigid lid configuration 



0.2 



0.4 



0.6 



0.8 



Figure 3. Dependence of the critical 
ratio S c on the density ratio 7, for both 
the rigid lid and free surface configura- 
tions. 
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Magnitude of the deformations. Depending of the parameters 7 and 8, we are able to compare the 

magnitudes of the respective amplitudes of the surface and the interface waves. Indeed, since 772 = S»=i u i e i.2 = 
Yl V±,±: an d Ci = Vi + V2 — J2i = i u i( e i,i + e ia) = Yl C±,±j one sees immediately that the surface and interface 
deformations, for each mode, are proportional, and satisfy 



V±,± 

C±,± 



e, 2 



- 1 



ei,i + e» 9 



In that way, one deduces that for the fast mode, the surface deformation is always bigger than the interface 
deformation, and the ratio tends to zero when 7 — > with S > 1, or when <5 — > 00. 

Meanwhile, as remarked in [22], the amplitude of the surface deformation is bigger than the one of the 
interface for the slow mode if < S < 2(1 — 27), and conversely if S > 2(1 — 27). Moreover, the ratio tends to 
zero when 7 — > with 6 < 1, and tends to 00 when 7 — > 1 or 5 — > 00. 

Thickness. In addition to forcing the polarity of the solitary wave, the ratio jM- is also related to the thickness, 
or the wavelength of the wave. Indeed, defining the thickness of a wave as in [34] by 



1 f + °° 



u{x)dx, 



one obtains for the function (32): l(u) = ^ = \ym- 

As an immediate result, and since we know the ratio between the magnitude of the deformations at the 
surface and the interface, the thickness of the deformations are identical at the surface and at the interface: 



A3 



1 



= 1. 



One can now compare the thickness of the different wave modes, and the ones in the rigid lid configuration, 
for waves of same heights. One computes in Figure 4 the ratio * %$-, ./tt and \/a? f° r 7 e (0, 1), and 
(5 = 1/2,1,2. 
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Figure 4. The different thicknesses, in the rigid lid and free surface configuration, for the 
parameters 7 € (0, 1) and (a) S = 1/2, (b) 8 = 1, (c) S = 2. 
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Of course, the thickness of the solitary wave tends to infinity when the depth ratio approaches its critical 
ratio 5 = 5 C , so that the waves predicted in the rigid lid and in the free surface configurations are excessively 
unlike. Additionally, one sees that when 5 is small, then the thickness of the internal waves predicted in the 
free surface configuration is largely different from both the slow mode and fast mode thicknesses, except in the 
limit cases 7 — > and 7 — > 1. 

When 7 <C 1 and 6 < 1, one has at the same time the similitude of the thickness and the velocities of the 
waves in the rigid lid configuration, and the fast mode waves in the free surface configuration. However, the 
amplitude of the surface deformation is bigger than the one of the interface for the fast mode waves, so that the 
rigid lid is not valid hypothesis. This resemblance has to do with the fact that both the fast mode waves and 
the waves in the rigid lid configuration converge to waves of the one-layer problem, when 7 — > and S — > 1~. 
Conversely, when 7 ~ 1 and when S 3> 1, the characteristics of the waves in the rigid lid configuration resemble 
the slow mode waves ones. We explore thereafter the validity of this hypothesis for different situations. 

3.4.2. Evolution of the rigid lid hypothesis 

Now, we restrict ourselves to initial data that are compatible for both the free surface and rigid lid configu- 
rations, and compare the different evolutions of the solutions. Following Section 3.3, if the initial data with the 
rigid lid assumption is (rj°,v°), then the corresponding initial data in the free surface case (in the limit e — >• 0) 
is U° = (—i] , r]° , —_ rgV°, — r^w ))- Therefore, one has initially «ii t=0 = e^ • SqU°, so 

'fe.fclt-o = Z, e ^| t= „ = e,, 2 I (1 - 7)^,377 + — — I , 

G\fc)=(±,±) i=i v 1 / 

0,fe| t=0 - Ki + ei, 2 ) I (1 - 7K 2 t?° + 

(j,fe)=(±,±) v ' 

Using the values given in Remark 3.5 page f3, this reads 

C° -C-, k ((1 -An I i ^ \±J*^1*A S\ 

Since the KdV equation preserves mass, knowing the size of the initial data allows to predict the significance of 
the waves. In particular, the rigid lid hypothesis will be valid for long times only if |QL ± I ^ \v± ± I- We give in 
Table 1 the different behavior of these variables, in the limits 7— > 1, 7 4 0, 5 ^ 00 and 5 — >• 0. As we can see, 
the rigid lid hypothesis will be valid for long times when 7 ~ 1, or when i > 1. For each of these cases, one 
sees that the main deformation comes from the slow mode waves, which correspond to the waves predicted by 
the models in the rigid lid configuration. 

As a specific example, when the initial data has zero velocities (that is to say v° =0), then we are able to 
compare straightforwardly the different magnitudes of the four waves. Indeed, one deduces from the previous 
calculations that when v° — 0, one has 



. fc(l-7) 4-1 _ fcO-7) 

J < fe 2S(c 2 + -c 2 _) c 2 k Vl ana ^ fc ~ 25{c\ -c 2 _) V ' 



,0 



Consequently, the four different waves have the same weight at the surface. The situation is more sophisticated 
at the interface, and one has eventually 

\v±,+ \li > \v±,-\li if <5 < 1-27, 
\v±,+ \l 2 < \v±-\l 2 if S > I-27, 
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S — > 00 
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\rP ± qv° 
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h(Tl°T^v°) 


Cl + 





1 ^0 1 _ 7 .o\ 


(t-,0 _i_ ,„,0\ 





1(1 -7)*?° 


2(^-1) ^ ±,w ^ 


2(1-5) ^ ±tU j 


Cl- 





1 


1 

1 „0 _|_ _,.o\ 





^((1-7)^ T^« ) 


2(5 - 1) ^ T W > 


2(1 - 5) t ?/ ± TV) 



with ^ 



H an di/ = 



7-j-l 2 % /(5+l)(l- 7 )' ' 7 ^0 *T """ ' VS(S-1) ' 

Table 1. Initial magnitudes of the different waves at the surface and at the interface, for an 
initial with a flat surface, in the limit cases. 



so that the fast mode wave is more significant than the slow mode wave when 5 < 1 — 27, and conversely 
otherwise. In the limits 5 — > 00 and 7 — > 1, the magnitudes of the fast mode waves tend to 0, so that the energy 
is only shared by the slow mode waves. Meanwhile, in the limit 5 — > 0, the significance of the fast mode and 
the slow modes respectively tend to -^-\t] \l2 and ^\rf\ L 2. Finally, in the limit 7 — > 0, then the magnitudes of 
the fast mode waves tend to when S > 1, and in the contrary carry all the energy when 6 < 1. 
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Figure 5. Magnitude of the deformations: the fast mode internal waves and the surface waves, 
when compared with the slow mode internal waves. Ratios for 7 <E (0, 1) and (a) 6 — 1/2, (b) 
S = 1, (c) S = 2. 

We plot in Figure 5 the different ratio of magnitudes, with the slow mode internal waves chosen as reference; 
that is to say 



IC±,+ I 



L 2 



Ki 



L 2 



and 



\V±,+ \l 2 



\V±-\l 2 \V±-\l 2 \V±-\l 2 

The rigid lid hypothesis is valid for small values of these ratios. Again, one sees that it occurs only when 7 ~ 1, 
or when S is big. This fact has already been addressed, for example in [25,34], but its precise confirmation has 
never been exposed, to our knowledge. For example, if 6 = 1/2 and 7 = 0.8, one has 






±.+ \L 2 



IC±,- 



L 2 



1 



|L 2 



\V±.-\l 2 



and 



\v±,+ \l* _j_ 
\v±,-\l 2 ~ 10 ' 



so that the rigid lid assumption is inaccurate. 
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3.4.3. Summary 

Let us summarize in Figure 6 the different results 
we obtained concerning the dependence of the behav- 
ior of the KdV approximation, depending on the pa- 
rameters of the problem. 

When 6 is above the plain curve (6 > 6 C ), then 
we know that the slow mode solitary waves will be 
of elevation type at the interface and of depression 
type at the surface, and conversely if 6 < 6 C . The fast 
mode solitary waves are always of elevation type at 
the surface and at the interface. 

Above the dashed line, the interface deformation 
is more important than the surface one for the slow 
mode solitary waves, and conversely if 6 < 2(1 — 2j). 
As for the fast mode waves, the surface elevation is 
always bigger than the interface elevation. 

The dash-dotted line concerns waves created by an 
initial data with zero velocities, and with a fiat sur- 
face. In that case, the fast mode waves will be smaller 
than the slow mode waves above the line S — 1 — 27, 
and conversely below. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

y 

Figure 6. Properties of the KdV ap- 
proximation, depending on the depths 
ratio 6 and the density ratio 7. 



We see that for big values of 7 and/or 6, the KdV approximation of the two-layer problem with a free surface 
gives a solution that resembles the interface problem with a rigid lid: the fast mode waves are smaller than the 
slow mode waves, and the magnitude of the deformation of the surface is of less importance than the deformation 
of the interface. On the contrary, when 7 and 6 are small, then the solutions of our problem, when considered 
at the interface, are comparable to solutions of the one-layer water wave models. 

The dots on Figure 6 represent the pair of parameters (7, 6) for which numerical simulations have been 
computed in Section 4.2. We have first computed the solutions of the symmetric Boussinesq/Boussinesq model 
and the KdV approximation, for both for the case of solitary waves, and zero velocities-flat surface initial data. 
The results are plotted for parameters corresponding to A (7 = 1/4, 6=1) in Figures 8 and 9, and to B 
(7 = 1/4, 6 = 2) in Figures 10 and 11. Then, we have compared the KdV approximation in the rigid lid and free 
surface configurations. Figure 12, 14, 13 and 15 corresponds respectively to the points A, B, C (7 = 9/10, 6=1) 
and D (7 = 5 = 1/4). 



4. Numerical comparison 



4.1. The numerical schemes 



This section is devoted to the numerical comparison between the different models displayed in this article, 
namely the symmetric Boussinesq/Boussinesq model (5), and the KdV equations (24). We first provide a numer- 
ical scheme for generic KdV equations that can easily be adjusted for the uncoupled KdV approximations (24) 
and (31), and its adaptation to the Boussinesq/Boussinesq system (5). 

Each time, we use a Crank Nicholson scheme and replace the costly numerical treatment of the nonlinear 
term by a predictive step. This method has been introduced by Besse and Bruneau in [4], justified in [3], and 
used in the water wave framework by Chazel in [10], and more recently by Durufle and Israwi in [17]. The 
method is formally of order two in space and time, which is confirmed by the simulations, and appears to be 
unconditionally stable. 



TITLE WILL BE SET BY THE PUBLISHER 27 

4.1.1. The KdV equation 

We present here the numerical scheme for the generic KdV equation: 

dtU + cd x u + Xud x u + nd x u. (33) 

First, we use the following semi-discretized in time equation: 

"" +1 d p" + cd x ( M " +1 2 +M " ) + A (au^d x (^^) + (1 - a)d x u"+^ (^^)) 
with a £ [0, 1], and where u n + 1 / 2 is a predictive term defined by 

u n+l/2 i u n-l/2 



(34) 



2 

The scheme takes advantage of the two possible discretizations of the nonlinear term ud x u, that is to say 
u ti+i/2q^ I u — +u \ an j u — +u q^ ^ u n+i/2^ ^y introducing a parameter a S [0, 1] and taking a convex com- 
bination of these possibilities. 

It is easy to check that, in order to preserve the semi-discrete L 2 -norm, one has to choose a = 2/3. As for the 
spatial discretization, we use the Crank-Nicholson scheme, adjusted so that the discrete L 2 -norm is preserved. 
This leads to the final discretization: 

u? +1 -u? . („ «" +1 + «" \ \(( <tl + u?.m + <+} + ul 



dl 



>„ + <r r sn ( Ci h=^ ) +M ( ft *^ . „, (35) 



with Di and D 3 the classical centered discretizations of the derivatives d x and d x with periodic boundary 
conditions. The scheme is given at each step by (35) and (34), with a simple explicit scheme for the first half 
step. One can then check (see [17, Theorem 2]): 



vneN, 5>n 2 = y, 



u°\ 2 . 



4.1.2. The Boussinesq/Boussinesq system 

We use the same ideas as in the previous section for the discretization of the Boussinesq/Boussinesq system. 
Even if the L 2 norm is not preserved by (5), and no simple quantity either, we decide to use the same parameter 
a = 2/3 as in the previous section. This leads to the following discretization of the spatial nonlinear term: 

^(u)d x u ~ \^{u^/ 2 )d x ( ^" +1 2 +t/ " ) + \^ ( ^" +1 2 +L/ " ) d x U^'\ 

As for the nonlinear term in time, we simply use 

' U n+1 + u r - 



Si(U)dtU ~ Si(f7" +1 / 2 )d t 
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We also need to construct the linear mappings Si with values on M^M), such that 

VU,VeR 4 , Zi{U)V = Si(V)J7. 

This finally leads to the following scheme: 

u? +1 -up ( U n+1 + U n \ ( U n+1 + U n ' 



(So + eS.iUr 1 ' 2 ) - eS 2 D 2 ) U < ~ U < + S ^i 2 + - - e (x 2 D 3 



/ rrn+l/2 ,-,.71+1/2- 

, l Ei ( K +1/2 + %+1 X 1 - 1 ! ( Or ^— ) W) 

Ei (lW" +1/2 




where Di, Z?2 and D3 are block-diagonal, with the classical centered discretizations of the derivatives d x , d 2 
and <9 3 (with periodic boundary conditions) as diagonal blocks. 

The scheme is given at each step by (36) and (34), with a simple explicit scheme for the first half step. 

4.1.3. Validation of the numerical method 

As said previously, the method is formally of order 2 in space and time, which is confirmed by the simulations, 
and appears to be unconditionally stable. In order to validate the schemes, we use the known solitary wave 
solutions of (33), expressed as follows: 

M 

u(t,x) = ) (37) 

cosh (k(x — cJtj) 

with c' = c+ ^, k = i/H^, and M arbitrary. 

Therefore, we are able to construct an initial data that will lead the solutions of the KdV approximation 
defined in Theorem 3.1 to be steady traveling waves. Indeed, if we set 7 € (0, 1) and 6 > 0, the coefficients of the 
KdV equations solved by the KdV approximation are explicit and given in Remark 3.5 page 13. Consequently, 
we choose Mi, . . . , M4 e R, and set U° = X)i=i u i( x i 0) e i; with Ui of the form given by (37), with the given 
parameters, and the uncoupled KdV approximation will be given by U(t,x) = Ylj—i Ui(et,x — Cji)ej. 

Unfortunately we cannot exhibit such an exact solution for the symmetric Boussincsq/Boussincsq system (5). 
In order to validate the scheme (36), we plug the solution of the uncoupled KdV approximation in (5), and 
obtain a forcing term F(t,x), so that U(t, x) is solution of the modified system 



(So + c(5i(J7) - S 2 dl))d t U + (S + e(2i(E0 - Z 2 d 2 x ))d x U 



F. 



It is trivial then to modify the scheme (36) by adding a forcing term F™ = F(t n ,Xi). 

We present in Table 2 the results that we obtain using the scheme (35) and (36), for several values of spatial 
and time discretization steps dx, dt, and different values of e, and for times T = 1/ e. The relative errors are 
computed in the discrete L 2 norm. These results allow to validate the schemes proposed. 



dx dt L T e KdV scheme Boussincsq scheme 
0.01 0.01 120 5 0^ 9.6317.10" 5 9.8719.10~ 5 

0.02 0.02 120 10 0.1 7.7094.10~ 4 7.9861. 10" 4 

0.05 0.05 120 20 0.05 9.5663. 10" 3 9.8587.1Q- 3 

Table 2. Numerical errors of the KdV and the Boussinesq schemes. 

In the following simulations, we always choose dx = dt = 0.01. 



TITLE WILL BE SET BY THE PUBLISHER 



29 



4.2. Numerical results 

We compute our schemes for two different different forms of initial data. The first leads to solitary waves as 
exact solutions of the KdV approximation (thus of the form (32)), and is therefore exponentially decreasing in 



space. The other initial data consist in deformations of the interface of the form 



M 



(therefore, they do not 



Vl + kx 2 

satisfy the spatial rapidly decreasing assumption for the better convergence rate in Proposition 3.f2), as the 
surface is flat, and initial velocities are zeros. All the forthcoming results are expressed in non-dimensionalizcd 
variables. 

We first look at the behavior of the relative difference between the solutions of the different models for 
fixed time, and for different values of e. In Table 3, we provide the difference between the original Boussi- 
nesq/Boussincsq model (4) and our symmetric Boussinesq/Boussinesq system (5), at time T = 1 and for differ- 
ent values of e. The same results are given for the comparison between the symmetric Boussinesq/Boussinesq 
system and the KdV approximation in Table 4. As predicted, when compared with the solution of the symmet- 
ric Boussinesq/Boussinesq model, the deviation of the original Boussinesq/Boussinesq system at fixed time is 
of order 0(e 2 ), whereas the error of the KdV approximation is of order 0(e). This supports our choice of the 
symmetric Boussinesq/Boussinesq system (5), as an equivalent model for system (4). 



dx dt 



L T 



0.01 0.01 120 I 0.1 
0.01 0.01 120 1 0.05 
0.01 0.01 120 1 0.01 



relative error of the solutions with the initial value: 
M „ n M 



U° 



cosh(/cx) 



U° = 



1.4859.10 



T 



y/1 + (kx) 



2.4785. 10" 



TT 



5.2714.10~ 4 
3.1656.10~ 5 



1.0882.10" 3 

7.1958. 10" 5 



Table 3. Relative error between the solutions of the symmetric Boussinesq/Boussinesq system 
and the original Boussinesq/Boussinesq model. 



dx dt 



L T 



0.01 0.01 120 1 0.1 
0.01 0.01 120 1 0.05 
0.01 0.01 120 1 0.01 



relative error of the solutions with the initial value: 
M „ n M 



U° = 



cosh(kx) 2 



U° = 



4.0805. 10" 3 
2.2785.10- 3 
5.2496.10~ 4 



y/1 + (kxf 



5.9295. 10"^" 

2.5506.10~ 3 

3.5434.10- 4 



Table 4. Relative error between the solution of the symmetric Boussinesq/Boussinesq system 
and the solution of the KdV approximation. 



Then, we give a numerical confirmation of the convergence rate obtained in Proposition 3.12. It is stated 
that for any initial data, the difference between the solution of the Boussinesq/Boussinesq system and the KdV 
approximation is bounded by 0(e) over times of order 0(l/e) if the initial data is sufficiently decreasing in 
space, and 0(ey/i) otherwise. 

We plot in Figure 7 the difference (in the discrete L 2 norm) of two solutions, obtained respectively by the KdV 
scheme and the Boussinesq scheme. The first plot concerns solitary waves of the form (32), thus exponentially 
decreasing in space, and the second is given by a zero velocities-flat surface initial data, with the interface of 
the form , — . 

^l + (kx) 2 

For both configurations, we set 5 = 1, 7 = 1/4, and simulate at the different values e = 0.1, 0.05, 0.025, 0.01, 
throughout time T = 1/e. One sees that the results match the theory, and that the decreasing in space at 
infinity is indeed of great concern. 
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0.045 

0.04 
0.035 

0.03 
0.025 

0.02 
0.015 

0.01 



•A 


— e=0.1 
---e=0.05 
E=0.025 

-— e-0.01 


I 


■ 




(a) solitary waves 



(b) zero velocities-flat surface initial value 



Figure 7. Relative error between the KdV approximation and the symmetric Boussi- 
nesq/Boussinesq system: (a) for solitary waves, exponentially decreasing in space (b) for zero 
velocities-flat surface initial data, slowly decreasing in space. 



The Figures 8 and 9 provide the snapshots of the previous simulations for e = 0.1, and at times t = 20 and 
t = 40. The results of the two schemes are plotted in the same figure, as well as the initial data and the (ten 
times) emphasized error at the interface, for readability. We only show the right side of the simulations, as the 
left side is obtained symmetrically. 

We see that in the case of a solitary wave, the Boussinesq/Boussinesq system produces an almost perfect 
soliton, that is predicted by the KdV approximation (Figure 8). The difference between the two solutions is 
located in space. However, when the initial data is not sufficiently decreasing at infinity, there is a big qualitative 
difference between the solution of the Boussinesq model, and the KdV approximation (Figure 9), since their 
difference is not exclusively located in the area of the solitary waves. Therefore, apart from a good initial 
localization in space, there is a significant interaction between the traveling waves of different wave modes, that 
is captured by the Boussinesq/Boussinesq system, and not by the KdV approximation. 

The same simulations are produced for different parameters, i.e. 5 — 2 and 7 = 1/4, and are displayed in 
Figures 10 and 11. We see that the exact same phenomenon appears. 

Finally, we present from Figures 12 to 15 numerical simulations of the KdV approximation, in both the 
free surface and rigid lid configurations. The simulation of the Boussinesq/Boussinesq models leads to almost 
identical results, so that we do not plot them for readability. Again, we set e = 0.1, and the different times of 
the snapshots are t — 20 and t — 40. In Figures 12 and 13, we choose the parameters 5=1, and respectively 
7 = 1/4 and 7 = 9/10. In Figures 14 and 15, we choose the parameters 7 = 1/4, and respectively 5 = 2 
and 5 = 1/4. Each time, the initial data consists in a bell curve in the interface, with no velocities and with 
a flat top. One clearly sees that, as discussed in Section 3.4, the rigid lid assumption is satisfactory as an 
approximation of the free surface, only in the case of small difference of densities (7 ~ 1), or when the depth 
ratio becomes large. 
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Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 




Initial data 
-Bouss/Bouss system 
- KdV approximation 
-errorxlO 



10 



20 



(a) t = 20. e = 0.1, 5 = 1,7 = 1/4. 



(b) t = 40. e = 0.1,(5 = 1,7= 1/4. 



FIGURE 8. Solitary wave solution of the KdV approximation, and the symmetric Boussi- 
nesq/Boussinesq model, at times (a) t=20, (b) t=40. 





Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 





" Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 



50 



50 



(a) t = 20. e = 0.1, 6 = 1,7 = 1/4. 



(b) t = 40. e= 0.1,(5 = 1,7 = 1/4. 



Figure 9. Solution of the KdV approximation, and symmetric Boussinesq/Boussinesq system 
for a zero velocities-flat surface initial value, at times (a) t=20, (b) t=40. 



32 



TITLE WILL BE SET BY THE PUBLISHER 



f 


1 


I A /v 


■ 


At, s\h 


r v 




11 v 




Initial data 

— Bouss/Bouss system 

— KdV approximation 

— errorxlO 



" Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 




(a) t = 20. e = 0.1, 6 = 2,7 = 1/4. 



(b) t = 40. e= 0.1,(5 = 2,7 = 1/4. 



Figure 10. Solitary wave solution of the KdV approximation, and symmetric Boussi- 
nesq/Boussinesq model, at times (a) t=20, (b) t=40. 




Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 



50 




" Initial data 
-Bouss/Bouss system 
-KdV approximation 
-errorxlO 



(a) t = 20. e = 0.1, S = 2, 7 = 1/4. 



(b) t = 40. € = 0.1,8 = 2,7 = 1/4. 



FIGURE 1 1 . Solution of the KdV approximation, and symmetric Boussinesq/Boussinesq system 
for a zero velocities-flat surface initial value, at times (a) t=20, (b) t=40. 
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Initial data 
-Free surface configuration 
- Rigid lid configuration 





t*\ 


■ 




- 




Initial data 

— Free surface configuration 

— Rigid lid configuration 



40 



50 



(a) t = 20. e = 0.1, S = 1, 7 = 1/4. 



(b) t = 40. e= 0.1,5 = 1,7= 1/4. 



Figure 12. Solution of the KdV approximations, in both rigid lid and free surface configura- 
tions, with two fluids of highly different densities, at times (a) t=20, (b) t=40. 





Initial data 
-Free surface configuration 
-Rigid lid configuration 





Initial data 
-Free surface configuration 
-Rigid lid configuration 



40 



50 



40 



50 



(a) t = 20. e = 0.1, S = 1,7 = 9/10. 



(b) i = 40. e = 0.1,<5= 1,7 = 9/10. 



Figure 13. Solution of the KdV approximations, in both rigid lid and free surface configura- 
tions, with two fluids of near equal densities, at times (a) t=20, (b) t=40. 
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Initial data 

Free surface configuration 

Rigid lid configuration 





Initial data 
-Free surface configuration 
-Rigid lid configuration 



40 



50 



40 



50 



(a) t = 20. e = 0.1, 5 = 2,7 = 1/4. 



(b) t = 40. e= 0.1,5 = 2,7 = 1/4. 



Figure 14. Solution of the KdV approximations, in both rigid lid and free surface configura- 
tions, with the upper fluid of greater depth, at times (a) t=20, (b) t=40. 



Initial data 
-Free surface configuration 
-Rigid lid configuration 



100 




Initial data 
-Free surface configuration 
-Rigid lid configuration 



100 120 



(a) t = 20. e = 0.1,8 = 1/4, 7 = 1/4. 



(b) t = 40. e = 0.1, 5 = 1/4, 7 = 1/4. 



Figure 15. Solution of the KdV approximations, in both rigid lid and free surface configura- 
tions, with the lower fluid of greater depth, at times (a) t=20, (b) t=40. 
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Appendix A. Proof of Proposition 2.6 

The proof is made of three steps. First, we introduce an energy of the system, and obtain an a priori estimate 
on this energy. Then, using this estimate and regularization operators, we prove the existence of a solution of 
our problem. Finally, using the energy estimate on the difference of two solutions, we get the uniqueness of the 
solution. 

A.l. Energy estimate 

In the following, we denote by U a solution of (5) on [0, T], with U G i°°([0, T}; H s+1 ) 4 . When we multiply (5) 
by A s , the system becomes 

(s + eSi(E0 - eS 2 dl)8 t K s U + e[A s , S^U^U + (s + e£i(tf) - eZ 2 %)d x A'U (3g) 

+e[A s ,S 1 (t/)]a rc C/ = 0. 

We then introduce the energy associated to the system 

Definition A.l. Let U <G H S+1 (M. 4 ). We define the energy of the function U associated to the system (5) as 

E S (U) = 1/2(S A S U, A S U) + e/2(S 1 (U)A s U, A S U) + e/2(S 2 A s d x U, A s d x U). 

If there is no risk of confusion, we simply write E s . 

From Assumption 2.3, So and S 2 are definite positive (with eigenvalues depending on 7 and S). Moreover, 
since s > 1/2, one has by Sobolev embeddings |«Si(/7)| Loo < C |C/| ff3 . Hence, there exists a = C ( n- v ^+|) 
such that if e £/ „. < ^ r \ , then 

-\U\ 2 H s+i < E S {U) < a\U\ 2 HS+1 . (39) 

Let us multiply (38) on the right by A S U and integrate. One obtains 

^ s -e/2(5 1 (a t C/)A s t/,A s C/)- e ([A s ,5 1 (C/)]9 t t/,A s C/) + e /2((S 1 (a rc C/)A s t/),A s t/) (4Q) 

-eflA'.EiCEO^l/.A'EO. 
Now, one has thanks to Cauchy-Schwarz inequality and Sobolev embeddings 

KSi^jyjA'^A"^)! < ^{dJJ^^K'U^ < C \U\ 3 H3 . (41) 

We then use classical Kato-Ponce commutator estimate [23]: 

Lemma A. 2. (Kato-Ponce) For s > 0, if f E H s and g e iJ s_1 7 then one has the estimate 

|[A S ,/]5| L2 < Co\f\ H s\g\ L oo +Co\d x f\ LOO \g\ Ha _ 1 . 
In our case, it leads to the commutator estimate 

|([A 8 ,£i(17)]0 x tf,A 8 17)| < laA'.Si^a^yA^I^Col^l^.. (42) 

In order to deal with the nonlinear terms with a time-derivative dt , we will need the following Lemma, using 
the elliptic form of the operator P e defined in (6): 

Lemma A. 3. Let P e (U, d) and Q e (U, d) the differential operators defined in (6) and s > 1/2. Then there exists 
Cq = Cq( ,i_ \ ,S+ -r) such that ifc\U\„ a < 1/Co, one has the following properties: 
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(1) P e (U,d) : H 1 — )• H^ 1 is one-to-one, and onto. 

(2) For k £ [0,s], j/one /ias P e (U,d)V £ /f*" 1 , tfien V £ H k+1 . 

(3) For 1/2 < fc < s and V £ if fe+1 , one /i«S 

1/Co\V\ h h +1 < (P € (U,d)V,V) < Co\V\ Hi+1 . 

(4) For 1/2 < k < s, the operator P € (U,d x )~ 1 Q e (U,d x ) is bounded H k+1 — > H k+1 , uniformly with respect 
to e. 

Proof. (1) Using the fact that So and S 2 are symmetric definite positive (with eigenvalues depending on 7 

and S) in the formulation of P e , it is obvious to see that one can choose Co = Cq{ ,^_ -, ,5 + 4) > 

such that for e|[/| LOO < 1/Co and any V £ F/ 1 , 

(P e (f/,d)V,F) = (S V,V) + e(S 1 (U)V,V) + e(S 2 d x V,d x V) > £-\V\ 2 H1 . 

Oo " 

In the same way, the bilinear form a(V, W) = (P e (U, d)V, W) is clearly continuous on H 1 x H 1 : 
{P e {U,d)V, V) = (SoV, V) + e(5i(C/)V, V) + e{S 2 d x V,d x V) <C \V\\,. 

Using Lax-Milgram lemma, we obtain that for all F £ H~ , there exists a unique V £ i? 1 such that for 
all W £ H 1 , a(V, W) = (F, JU), and hence there exists a unique variational solution of 

P e (U,d)V = F. 

(2) We will prove the second point by induction. The result is known for k = from the previous point. 
Then, we remark that if P e (U, d)V = W with W £ L 2 (and U £ ff 1 ), then 

eS^V = S V + eS 1 (U)V - W. 

Therefore, since S 2 is invcrtible, we get <9 2 U £ L 2 , and thus V £ H 2 . The second point is therefore 
proved for k = 1, and the intermediate values fc £ [0, 1] follow by interpolation. 

Let now assume that the result is known for a k £ N* and let P e (U, d)V £ H k C H k ~ 1 . From the 
induction hypothesis, one has V £ H k+1 , and moreover 

eS 2 d 2 A k V = S A k V + eA k (S 1 (U)V) - A k P e (U,d)V. 

We know that A k V £ F 2 and A k P e (U, d)V £ F 2 . Moreover, since U £ H s and s > k > 1/2, we know 
from the Sobolev embedding that A k (S 1 (U)V) £ L 2 and |A fc (Si(17)F)| 2 < C |f7| Hfc \V\ Hk . Therefore, 
d 2 A k V £ i 2 , and as we know from the induction hypothesis that V £ H k+1 , one has indeed V £ H k+2 . 
The second point is thus proved for all k £ N(~l [0, s], and we obtain the result for < k < s in the same 
way, starting the induction from k — [fcj £ [0, 1). 

(3) Since So and S 2 are symmetric definite positive, one has for t\U\ „ 3 < 1/Co with Co = Co( ^3 y , o~ + i), 
and for any V £ H k+1 , 

(A k (P f XU,d)V),A k V) = (S A k V,A k V) + e(A k (S 1 (U)V),A k V)+e(S 2 d x A k V,d x A k V) 

1 I i2 1 I i2 

— 7^ Lff fc T^ 6 H-fc+i ■ 
Go Oo 

The second inequality is straightforward. 
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(4) Let W G H k+1 . It is obvious that, since s > k > 1/2, one has 

Q e (U, d x )W = (S + cSi(17)) W - eX 2 d 2 x W e ff*" 1 . 
Hence, we know from the previous points that V = P e (E/, d x )~ 1 Q e (U, d x )W € H k+1 , and 

i|^li*+i < (A fe (P £ (C/,a)F),A fc T/) < C \V\ 2 Hk+1 . 

Therefore, one has: 

\V\ 2 Hi+ i <C (A k (Q t (U,d)W),A k V) 

<C (A k (^o + e^i(U))W),A k V)+C (d x A k J: 2 W),d x A k V) 

< C (7, J, |^| ff *)(|W| fffc |y| fffc + e|W| fffc+1 |V| fffc+1 ). 

Finally, one has 

M^ +1 <l(C \w\ Hk - \v\ Hk f + i(C |T^| fffc+1 - \v\ Hk+1 f + ^|^ fc+1 + \CI\W\ 2 H , +1 

<\\v\ 2 Hi+1 + \c%\w\ 2 Hi+1 . 

The operator P e (C7, d x )' 1 Q e (U, d x ) is thus bounded H k+1 -)• P e fc+1 by C5 2 , which ends the proof of the 
Lemma. 

□ 

Using this Lemma, we immediately obtain that for U <G H s+l satisfying (6) and the hypothesis of the Lemma, 
one has d t U € H s , and 

|s t [/| H| = Ip^cy-^^cy^t/l^ < c o |0xi7| H . < c |t/| ff|+ i. (43) 

Therefore, we can use the same calculations as in (41) and (42), and obtain 

|(^(a t C/)A s t/,A s C/)| < |Si(cW)| Loo |A s [/|' 2 < C \U\ 2 HS \U\ H;+1 (44) 

(since s — 1 > 1/2 such that |<9t£/| Loo < Co\dtU\„ s _ 1 < Co\u\ „, + i), and with Kato-Ponce theorem, 

|([A»,Si(EO]fttf,A'EOI < |([A s ,^i(C/)]9 t C/| L2 |A s C/| L2 

< ColC^I^-. l^e^lfi-— i ^ c 'o|^|^|l7| ff . + i. (45) 



Finally, one deduces from (39), (40) and the estimates (41)-(45) 

f t E s <eC \U\ 2 HS \U\ m 



4-E s < eC \U\ 2 „.\U\„. +1 < eC E 3 J 2 , 



with Co = Co( (Yz \ ,8 + |), and providing the fact that e|t/| Loofl - s < 1/Co- 
From Gronwall-Bihari's inequality and (39), it follows 

\U\ H!+1 < C E^ < Co f 2] ^ /2 < Co - ^IS f ( 46 ) 

l-C etE s > | t=0 l-O e|f \ H s+it 
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Thus, one sees that there exists Co such that if e\Uo\ H e+i < 1/Co, then one can choose T — T(Co) such that the 
smallness assumption e \U „ s < 1/Co remains valid for any t G [0,T/e]. Hence, in the frame of the proposition, 
one has the following estimate: 

\u°\ 

\ U \L-ao,T, e]i Hi») ^ C \-C e\U°\ H;+1 t' (4?) 



with Co = Co( T^z ) ,$ + j) an d T > 0, independent of e 

Let us note that t 
with the same bound. 

A. 2. Existence of a solution 



Let us note that by Lemma A. 3, we know that we can also control the time derivative \dtU\ Toon „ „, , „.-, 



We can deduce from the energy estimate (47) the existence of a maximal solution U € L t '°°([0,T/e]; H s+1 ) 
of (5) for any initial data U° <G H s+1 . We follow the classical Friedrichs proof, using the regularization operators 
defined thanks to the Fourier transform as below: 

V«eL 2 , V£eR, j>(£) = <p(v$v(S), (48) 

with (p a smooth numeric function with compact support, such that (f(0) = 1. These operators have the 
following classical properties 

Lemma A. 4. (1) J„ is bounded H s — > H s : there exists Co(s, v) such that for any v € H s , 

\J v v\ Ha <C {s,v)\v\ HS (49) 

(2) J u commutes with K s , and is a self-adjoint operator 

(3) There exists Cq independent of v such that 

|J^| 2 < C |v| 2 (50) 

We obtain then a solution of (6) as the limit of U v the solutions of 

d t u v + JvP^JvU^d^Q^JvU^d^j^Uv = o, (51) 

with U u]t=0 = U°. 

From Lemma A. 4, one has that (51) is an ordinary differential equation on the Banach space H s . Thus, 
thanks to Cauchy-Lipschitz Theorem, we know that there exists a unique maximal solution U v € C([0, T v ), H s ). 
Using the calculations of Section A.l and the properties of J u , one obtains the energy estimate 

1/7° I 



H <-~"l-C e\U°\ H:+1 t' 

with Co independent of v. Moreover, U v satisfies (43): 

\dtU v \ Hs _ t < C \d x U\ Hs _ t < C'o\U\ H3 . (52) 

Thus one can find T > independent of v such that the solution U v does not blow up for t e [0, T/e]. In 
particular, one has T v > T/e > 0. 
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Now, since (U v ) v> o is uniformly bounded in _L t ,oo ([0,T/e],_fff +1 ), one can extract a subsequence that weakly 
converges towards a function U € L t ' co ([0,T/e], iJ e s+1 ). We want to use Ascoli theorem, but we need the 
injection H S (M) C .£P _1 (R) to be compact, which is not true since R is unbounded. However, one can easily 
circumvent this problem, using weighted Sobolev spaces for example. Finally, one obtains that there exists 
U € C°([0,T/e),i7 s_1 ) such that U u converges strongly towards U as the subsequence v tends to 0. Also, by 
interpolation inequalities, U v converges strongly towards U in C°([0,T/e), H s ~ a ), for < a < 1. Then, since 
one can find a such that H s ~ a injects continuously in C 1 (R), one proves that Q e (J^U I/ ,d x )J^d x U v converges 
to Q e (U,d x )d x lI and J^P e (J u U„,d x )d t U converges to P e {U Vl d x )d t U as v — > 0. Hence, U is indeed a solution 
of (6). 

From Section A.l, we know that U € L t '°°([0,T/e]; _ff e s+1 ), and one can prove (see [40] XVI.1.4 for example) 
that Ue C a ([0,T/e);H s + 1 )f]C 1 ([0,T/e);H s ). 

A. 3. Uniqueness of the solution 

Let U 1 ,U 2 e C°([0,T/e);H s+1 ) n ^([O.T/e);^ 8 ) be two solutions of the Cauchy problem (5) with initial 
data Uu t=0 — U2\ t=0 — U . One can immediately check that R = U\ — U 2 satisfies 

(s* + eSi(E/i) - eS 2 d^d t A s R + (e + eEi(J7i) - eZ 2 d^d x A s R 

+e[A s ,S 1 (U 1 )]d t R + e[A s , Ei(E7i)]9 x iJ - eJ F, (53) 



with F = —A" I S\(R)dtU2 + Ei(_R)9a;f72 ) • Then, we can carry out the same calculations as in Section A.l on 
R, and obtain the equivalent energy estimate 

j f Es{R) < eC (|J7i| ff . + \U 2 \ H .)E„ 

with C = 00(^^,6+ i, |£/°| H|+1 ). 

From Gronwall-Bihari's inequality and the estimate (47) on U\ and U 2 , and since F s (i?)| f=0 = 0, one has 
immediately E S (R) =0 on [0,T/e], and finally U\ — U 2 . 
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